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Chapter 1 
The Electrostatic Field et 


1.1 Field Intensity An electrostatic field is described by specifying the 
electric intensity vector E at every point. £ is the force per unit charge. 
Thus, if a point charge q is placed at a point where the electric intensity 
is E, it experiences a force F given by 


F = 4E. Giny 


The intensity at P due to a single point charge q at O acts along OP — 
and has magnitude which is inversely proportional to OP*. Thus, if RS. 
OP = rand P is the unit vector along OP, then = 


GC RelA) 12)! 
cali =a pe ae pm sc)" 
In this equation, all quantities are measured in SI units (see p. 95) and ; 
b) = 885% 10~ 12 is called the permittivity of free space. 
It follows from equations (1.1) and (1.2) that the force exerted upon a 


charge q/ at P by the charge q at O is given by 


This is Coulomb's Law. re: 

The electric intensity due to any given distribution of charge can be 
calculated by vector addition of the intensities due to the individual 
charges. 


Problem 1.1 The surface of a disc of ebonite is uniformly charged. 
Calculate the electric intensity at points on the disc’s axis. 


Solution. Let o be the charge per unit area on the disc. Consider the 
contribution to the intensity at P (Fig, 1.1) of the elemental ring of radius F 
rand thickness dr. Any charge q at a point Q on this element generates an 
intensity g/4ne) PQ? at the point P in the direction QP. It is clear from 
the symmetry that the resultant intensity at P is directed along OP; 
this implies that only the component along OP of q/4n& PQ?, namely 
q.cos 6/4m@) PQ?, contributes to this resultant. Thus, the net contribution 
of the elemental ring is (5, q)cos 0/4me, PQ?, where 5) q is the total charge 
on the ring. But 2nr dr is the area of the ring, so that ) q = 2zor dr. 
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EM—B 


al o 
E=— i =—(l- 
mal sin 8 dO 2! Cos a), 


e a is the angle subtended by a radius of the disc at P. 


roblem 1.2 A very long straight rod is uniformly char; Calculate 
e electric intensity at a distance p from the rod. am 


Solution. The intensity at a point P (Fig. 1.2) must be directed along the 

cular OP to the rod. Let Q be a point on the rod distant x from 

and let QQ’ = dx be an element of the rod. Then, if q is the charge 

: length, the intensity at P due to the element has magnitude 

dx/4ne) PQ’ and is directed along QP. Only the component in the 

direction OP contributes to the resultant intensity at P. It follows by 
integration that this resultant is given by 


Since x = ptan0, dx = psec0 d0, PQ = psec, this integral is equ 


q {* 
E= cos 6 dO 
Ane, p )- in 


lent to 


Problem 1.3 An ebonite rod is uniformly charged, the total 
being q. It is placed along the axis of a disc of ebonite, also unifo 
charged with charge Q, so that one end of the rod is at the centre of 
disc. If the radius of the disc and the length of the rod are both a, show 
the rod is repelled by a force Qq(2—./2)/(2néq a”). 


Solution. Consider an element dx of the rod, distant x from the centr 
of the disc. The charge on the element is qdx/a; it lies in a field of in ‘“ 
Q(1—cos 0)/(2ne, a?) generated by the disc (see equation (1.4)), where 
is the angle subtended by a radius of the disc at the element. Thus, 
element is repelled along the disc’s axis by a force Qq(1—cos 0) dx/(2né a 4: \ 
Since cos@ = x/,/(a? +x), by integration with respect to x fromx=0 | 
to x = a, the resultant force on the rod is now easily found to be that 
stated, oe 


1.2 Potential When a point charge is taken from one point to anot 

in an electrostatic field, the work done by the field force acting on the 
charge is always independent of the path joining the points. Any field 
for which this principle holds is said to be conservative. 


2. it at 
aU ig 


z ‘ ane 3 
done by the field force when a point charge is taken from a 
‘toa fixed datum point D is called the potential energy of the ‘ 
‘The potential energy of unit charge at P is called the potential of the 
at P and will be denoted by V. 
is usual to choose the datum point to be at a great distance from the 
charges responsible for the field. In the case of the field due to a single 
ogg 4, the force exerted on unit charge distant r from q is q(4neyr); 
the work done by this force when the unit charge recedes to an 


~ 
& e distance is given by 


Eee 
V= herr (1.6) 
is the potential at distance r from q. It follows from equation (1.2) 
in this case, 
E = —gradV. (1.7) 
every electrostatic field is generated by point charges (electrons, 
protons, etc.), this equation is valid for all such fields. 


Problem 1.4 Calculate the potential of an infinite line charge q per unit 


th, 


_ Solution. The electric intensity has been calculated in Problem 1.2. 
is Perpendicular to the line charge and its magnitude is given by equa- 
1.5), Clearly, V depends only on the distance p from the line charge. 
by equation (1.7), the derivative of V in the direction of p must equal 

= —q/2néo p. Hence integrating, 


q 
v 75 log p. (1.8) 


The constant of integration is taken to be zero, implying that the datum 
p = 1. Since V-+ —co as p — on, the datum cannot be taken at 


infinity in this case. o 


__ A dipole comprises point charges +q and —q at points A, A’ (Fig. 1.3). 
Their distance apart A’A = d is assumed to be small and q is assumed to 
be correspondingly large so that the product qd is finite; this product 
he ‘is termed the moment of the dipole and will be denoted by m. The moment 
is usually taken to be a vector quantity, its direction being that of the 
_ dipole axis A’A; it will then be denoted by m. By summing the contribu- 


--eo) 


AA. 


Fig, 1.3 


Problem1.5 Calculate the field intensity due to a dipole. 


Solution. Let E,, Ey be the components of £ in the radial and transverse _ 
directions at P (Fig. 1.3) ie. along OP and perpendicular to OP. Then, by — 
equation (1.7), y l 

a __ lav ‘ 
B= =F Eom 73H" ( ‘: 
Substituting for V from equation (1.9). we find 
1 2mcosé 1__msin@ 
Ta a 4n, 

This result can also be expressed in vector form thus: We note that the f 
vector m/4néor° has components (m cos 0/4zéq r°, —m sin 0/4769 1°) in) 18 
the radial and transverse directions. It follows from equations (1.11) that “" iu 
the vector E+m/4ne)r*? has zero transverse component and radial 3 
component 3m cos 6/4ne, r?. Hence we can write 


vie. 


E, Ses (itt) 


loom 1 3mcos0 
 anay Panay 
a a) mer 
ata t 
Thus 
1 (3m.r om 
E= ae ( 5). (1.12) 


Problem1.6 A particle of mass M and charge q is free to move in a plane 
containing the axis of the dipole shown in Fig. 1.3. At t = 0, it is projected 
from the point r = a, 0 = a, with transverse velocity V and zero radial 
velocity. Prove that, at any later time t,r? = a? +kt?, where k is a constant. 
Hence find the condition for the particle to describe a circle. 


Solution. The energy equation for the particle is 


4 U) qm cos & 
M(?? +7262) 4 EOS = py? 4 ESE 

aM +70) 4neg r* 3 = 4néy a? 
Using the first of equations (1.11), the radial component of the particle’s 
equation of motion can be written 


M@r—r62) = md cos 
4neyr 


Eliminating 4? between these equations, we derive the equation 


+P? = 2k (1.13) 
where qm cosa 
4n@) Ma?” 

The left-hand member of equation (1.13) is d(r#)/dt. Hence, integrating 
under the initial condition * = 0 at t = 0, we find that r* = 2kt. Since 
r? = d(}r*)/dt, a further integration under the initial condition r = a at 
t = 0 yields the stated equation. 

The condition for a circular path is clearly k = 0, i.e. 

m = —2n& MV7a?/(q cos «). oO 


k =4V74 


Problem 1.7 Charge is uniformly distributed with density p over the 
interior of a hemisphere of radius b. Calculate the potential at a point on 
the axis of the hemisphere, on the side remote from the plane face, at 
distance y from the centre. If a particle of mass m and charge —2npb°/3 is 
released from rest at a point on the axis where y has a large value, find the 
velocity with which it arrives at the curved surface. 


a Be eee eee ae te a eT N 


Solution. If a thin disc of radius a has a uniform charge o per unit area, 
it follows from equation (1.4) that the electric intensity at a point on its 
axis distant 2(> 0) from the centre is given by 


pei eel ee | 
ar (2? +a7)* |" 
If V is the potential, then equation (1.7) requires that E = —dV/dz. Hence 
V= Z[erar-|. (1.14) 
(N.B. This makes V + Oasz—> 0) 


Dividing the hemisphere into disc-shaped elements of thickness dx as 
indicated in Fig, 1.4, x being the distance of an element from the centre O, 


Fig. 1.4 


it now follows that the contribution of such an element to the potential 
at P(OP = ») is 
dx 
oe {(y—xP + -?}8+x-y] = & [02 +b? —2yx)# +x]. 
by ‘0 if 
Integrating with respect to x over the range 0 < x < b, this gives a net 
potential at P of 
Le 24 ya\t— 23 — 3b? 3 
= = — [2b —2y? —3b?y+2b9]. 
ororarl te eam 


ing i i i i vo 

By expanding in ascending powers of 1/y, it may be verified that ; 
as y+ oo, At A, y = b, and V = pb?(4./2—3)/12e. Hence the loss in 
potential energy of the charge —2npb*/3 as it moves from a great distance 
to A is 2p2b5(4,/2—3)/18e9. This must be compensated by an equal gain 


7 ee - 
" Cer ee 


follows that : 
4 2 _ mp7b* ba 
= Fay 2-3. a 
Problem 1.8 An anode of a cathode ray tube is a plate bounded by 
concentric circles of radii a and Sa. It is to be assumed uniformly charged 
_ and non-conducting. An electron is stationary at a great distance from the 
Bee n its axis. It moves under the attraction of the plate, passing through 
centre with velocity v. Calculate its velocity v when at a distance 2a,/6 
from the plate, 


Solution. It follows from equation (1.14) that the potential due to the 
_ anode at a point on its axis distant z from its centre is given by 


V =< [le +250)— Je? +a%)]. 
; 


Thus, at z = 2a,/6 and z = 0, the electric potentials are ao/é, and 
2aG/e respectively. If —e is the charge on the electron, its potential energies 
at these two points are accordingly —aeo/e and —2aea/e. Since the 
fa ae total energy of the electron is zero, conservation of energy implies 
i at 

: Peart OS put 2 9; 

& & 


v 


m being the electron’s mass. It follows immediately that v = u//2. 0 


‘1.3 Gauss’s Law If S is a closed surface described in an electrostatic 

fi field and E, is the component of E in the direction of the outward normal 
at any point on S, the surface integral of E, over S is termed the flux of E 
out of S. Denoting the flux by N, we have 


af N=| Eds =| E.nds, (1.15) 
s s 


where n is the unit outwardly directed normal to S. 
J Gauss’s Law follows from the inverse square law (L. Marder, Vector 
PY Fields, p. 26) (equation (1.2)) and states that 

, N = Qty (1.16) 
% where Q is the total charge enclosed by S. 

By applying this result to a sphere of radius r which is concentric with 
a uniformly charged spherical surface, it is easily shown that the field 
intensity of such a charge distribution is given by 


E-7,. r>a; E=0, r<a, 


- of kinetic energy 4mv?, where v is the velocity of arrival at A. It now 


ye we , _ oo oui & ae 
att" ke a aro oF 
Here, is the total charge and £ is clearly in the radial direction. Thus, 
cere eed Gal onist wiki Gat Wilch ‘wou bs geome ME r 
whole charge Q were placed at the centre, : ; 
The corresponding potential function outside the sphere r = a is 


ed 
4neyr 


Vv 


(see equation (1.6). Inside this sphere, V must be constant and equal to — 
its value at the surface, viz Q/(47é, a). i 


Problem1.9 Charge is distributed uniformly throughout the interior of 
a sphere of radius a. Show that the magnitude of the field intensity within — 
the sphere is proportional to the distance from the centre. 


Solution. Considerations of symmetry indicate that E is radial. Take a i. 
Gauss surface in the form of a concentric sphere of radius r(< a). IfQisthe 
total charge, the charge contained by the Gauss sphere is Qr'/a*. The 
magnitude of E must be constant over this sphere and its flux out of the P 
sphere is accordingly 4zr?E. Gauss’s law now yields 4xr?E = Or Jeg a. 
Thus E = Or/(4ne, a°). aid 


Problem 1.10 A distribution of electric charge lies within an infinite , - 
circular cylinder of radius a. The charge density at a distance nA< a) from 
the axis is given by p = 3Q(a—r)/za*. Calculate the potential function. 


Solution. It follows from the axial symmetry that the intensity at any — 
point P distant r from the axis of the cylinder will be radial (Fig. 1.5). 
Construct a coaxial cylinder of unit length with its curved surface contain- 
ing P. If r < a, the net charge contained in this cylinder is 


Y yr? 
fi 2npr dr = 4 ra—r) dr = QF" (a—2r), 
° @ Jo a 


The flux of E out of the circular ends of the cylinder is zero and the flux 
out of the curved surface is 2nrE. Hence, by Gauss’s law, 
2nrE = Qr*(3a—2r)/éo a°). This gives E = Qr(3a—2r)/(2né a’). : 
If r > a, the total charge contained in the Gauss cylinder is Q and, — 
hence, E = Q/(2méyr). : 
If V is the potential, E = —dV/dr. Integration accordingly leads to the 
Tesults 
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Ri bat oat 4a). rea 


12m8 a? 
al — 3S otra) r>a. 
(Note: The datum point has been taken on the surface r = a, so that V 


vanishes on this surface.) o 


Fig. 1.5 


es The molecules of a conducting material contain charges which are free 
to move from molecule to molecule under the influence of any applied 

yi electric field. Thus, when a solid conductor is placed in an electric field, 
. electric currents will commence to flow in its interior and will only cease 
when the electric intensity over the interior has been reduced to zero by a 
redistribution of charge (the induced charge). Under static conditions, 
therefore, E = 0 within a conductor; a direct consequence of this is that 
V is constant over the conductor. In particular, the earth is an immense 
conductor whose potential is usually taken to be the datum with respect 
to which all other potentials are measured. Thus, the earth’s potential is 
zero and when a conductor is ‘earthed’ its potential is reduced to zero. 
Another consequence is that the interior of a conductor is uncharged; 
this follows since the flux out of any closed surface drawn inside the con- 
ductor is zero, and hence, by Gauss’s law, it can never contain any charge. 


10 


1.4 Laplace’s Equation Let S'be any closed surface drawn in an electro- 
static field and let I’ be the region enclosed by S. If p is the charge density, 
the total charge inside S is given by the volume integral | pdv, Hence, by 


Gauss’s law, 


[ B, ds = 05" pav. 
s e 
But, by the divergence theorem (L. Marder, Vector Fields, p. 43), 
| z,as = | divedo. 
Ss r 


It now follows that 


J (div E~pjeo) dv = 0. af 
r am 
Since this is true for every region I, the integrand must vanish. Thus, 
divE = p/eo. (1.47) 


This is the second fundamental equation of the electrostatic field. 
Eliminating E between equation (1.7) and (1.17), we obtain Poisson's — 


equation for the potential, namely, y 
div grad V = V?V = —p/éo. (08 
Thus, at all points where there is no charge, V satisfies Laplace’s equation, 


vv =0. (1.19) 
Any solution of this equation is called a harmonic function. 

It is assumed in equation (1.17) that p is finite and continuous. If the 
charge is distributed over a surface S asa film having infinitesimal thickness 
(e.g. S may be the surface of a conductor), this will not be the case. In these 
circumstances, let the region on one side of S be called region 1 and that — 
on the other side be region 2 and let m be the unit normal to S in the sense 
from 1 to 2 (Fig. 1.6). Then, if E,,, E2, are the components of £ in the | 
direction of n on the two sides of S, the flux of £ out of the small cylinder 
shown in the figure is (E,,—E,,)o, where is the area of its right-section 


/ 


Fig, 1.6 


h of the Gollahes odds acy nd w de oct 1 
the charge enclosed by the cylinder is cw. Thus, by Gauss’s law, 
Ey E yy = o/b. (1.20) 
particular, if S is the surface of a conductor and region 1 is its interior, 
E, = 0 and equation (1.20) reduces to 
E, = 0/€. (1.21) 
e subscript 2 can now be dropped.) 
It is further assumed that the potential V is continuous across a surface 
inuity, This is equivalent to the statement that any component 
gential to the surface is continuous. Thus, if the surface is that of a 
nductor, any tangential component of E just outside the conductor is 
and £ is accordingly normal to the surface everywhere; this also 
vs from the fact that such a surface is an equipotential. 
"For use in many of the later problems, it will be convenient to list a 
iber of well-known harmonic functions. 
If (p, $, 2) are cylindrical polar coordinates (Fig. 1.7) the following are 
ions of V?V = 0; 
V = Alogp+B, (1.22) 
= (Ap"+ Bp~"(C cos n@+D sin ng), (1.23) 
A, B,C, “sl are arbitrary constants and n is an integer. Since these 
solutions are independent of z, they give the same field in every plane 
el to Oxy, 
+ 9, $) are spherical polar coordinates (Fig. 1.7), the following is a 
lution of Laplace’s equation 


z 


Ve (arr) Pleo, 
where n is zero or a positive integer and P,(1) is Legendre’s polyn 
degree n; in particular 
Pou) = 1, 
Py) = 
Pu) = 430? —1), 
Ps(u) = Sp? 3p). 

More complex solutions can be derived from equations (1.22)-(1.24) if A 
addition (Principle of Superposition). we 
Problem 1.11 A solid conducting sphere of radius a is surrounded by 
concentric conducting spherical shell of internal radius b, which is earthed, 

If the inner sphere is raised to a potential ¢, calculate its total charge. 


Solution. Let P be any point in the field between the spheres and 
OP = r (Fig, 1.8). In this region, V must satisfy Laplace’s equation and be 


spherically symmetric about O. With n = 0, the solution (1.24) is indepen- 
dent of @ and, hence, spherically symmetric. Thus, we take 


V= ree 
r 


eis create tn becatey cat Y= = gover r =a, (i) V=0 
aes se ie 


Ate = 6, A+2 =0. 


_ Solving for A, B and substituting in equation (1.26), we find 


ap (b 
Clearly, the electric intensity is directed radially and its magnitude is 
accordingly given by 


Ve 


ov 1 
Thus, the flux of E across a concentric sphere of radius r is 
vam 4nabp 
pr b-a’ 
__ By Gauss’s law, this equals Q/ep, where Q is the charge on the insulated 
sphere. Hence 
ot _ Ane q abp 
i Q= ae eae (1.29) 


« Alternatively, the surface density of the charge on the sphere r = a can 
be calculated from equation (1.21) to be 

bb 

a(b—a) 

and this then yields Q as before. o 

i Problem 1.12 An earthed conducting sphere of radius a is,placed in a 


uniform field of intensity E, . Calculate the charge induced ona jisphere 
whose axis is in the direction of the field. Ps 


7 6 =6£E,n,= (1.30) 


Solution. Taking the centre of the sphere as origin O and the z-axis in 
the direction of the uniform field E, (Fig. 1.9), let r, @ be spherical polar 
_ coordinates of any point P of the field (the third coordinate ¢ will not be 
needed, since the field has axial symmetry about Oz). In the absence of 
the sphere, the uniform field has potential V satisfying the equations 
OV /éx = aV/éay = 0, AV/éz = —E,: thus 
4 V = —E,z = —Egrcos@ = —EprP,(cos@). 
, When the sphere is introduced, the charge induced on its surface causes 
a distortion of the uniform field, but this distortion will vanish at great 


14 
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Fig, 1.9 7 


distances from the sphere. Referring to equation (1.24), taking n = 1, we 
note that BP, (cos @)/r? is a solution of Laplace's equation vanishing at 
infinity. Assuming this to be the potential due to the induced charge, the 
resultant potential of the field becomes 


V= ~Ere0s0+%,cos0. (1.31) — 


It remains to satisfy the boundary condition over the sphere by appro- 
priate choice of B. This condition is V = 0 over r = a and requires that 
B = E,a’. Hence , 


a ; 
V=E, S-1)rcose 1.3) 


It has been shown that this potential function satisfies all the given 
conditions of the problem, but it is by no means self-evident that other _ 
potentials do not exist which also satisfy all these conditions. However, a 
uniqueness theorem can be proved (L. Marder, Vector Fields, p. 55) which 
asserts that this is not the case and equation (1.32) is accordingly estab- ik 
lished as the proper solution to our problem. Such uniqueness theorems Py 
will be very frequently appealed to in similar circumstances throughout 
the remainder of this book. On the surface of a conductor the conditions 
to be satisfied for uniqueness are that either the potential must take a 
given value or the total charge must be as given; if the charge is given, a 
further condition is that the potential must be constant (unknown) over 
the conductor. 

The radial component of E can now be calculated from equation (1.32), 
thus, 


Beatie (25 +1 )eo80 (1.33) 
15 


Putting r = a, this gives the normal component of E over the sphere to be 
3£E, cos 8. Equation (1.21) then shows that the surface density of the induced 
charge is ¢ = 329 Ecos, ie. the induced charge is positive over the 
hemisphere 0 < @ < 4m and is negative over the other hemisphere. By 
dividing the positively charged hemisphere into elemental rings by planes 
perpendicular to Oz, the total charge on this hemisphere is proved to be 

2na [* asin 0 dO = 3nae Eo. oO (1.34) 
Problem1.13 A dipole of moment m is placed at the centre of a spherical 
cavity of radius a in an earthed conductor. Calculate the field. 


Solution. Introducing spherical polar coordinates r, @ with pole O at 
the centre of the cavity and z-axis along the axis of the dipole (Fig. 1.10), 


Fig. 1.10 


we know from equation (1.9) that the potential due to the dipole alone is 
mcos 0/4néy r?. To this potential function must be added a term represent- 
ing the contribution of the charge induced on the wall of the cavity. This 
term must have axial symmetry about Oz and must be finite within the 
cavity; we shall assume it to take the form ArP,(cos@) = Arcos@ (the 
case B = 0,n = 1, of equation (1.24). Thus, the total potential within the 
cavity is given by are 
m 


Ve 
4ne, r* 


+ Arcos 6. (1.35) 
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It remains to satisfy the boundary condition V =0 over r =a. 
Equation (1.35) shows that V vanishes for all values of @ when r = a, 
provided A = —m/4ne, a°. Thus, 


m (@r 
vega (5-5) cos. (136) 
The radial component of the field is now calculated to be 
ov m 2a* 
E,= -7= ae (+1) cos 0. (1,37) 


Putting r = a, we find that the intensity at the wall of the cavity is 
—3m cos 0/47é, a°. Thus, the surface density of the charge induced on the 
wall is given by ¢ = —3mcos6/4za°. o 


1.5 Force on the Surface of a Conductor The film of charge on the surface 
of a conductor is situated in a field whose intensity decreases very rapidly 
across the film from o/é) just outside the conductor to 0 inside the con- 
ductor. The average field is accordingly o/2e, and the force per unit area 
applied to the conductor surface can be proved to be 07/29 directed along 
the outwards normal. 


Problem 1.14 A large plane conducting plate of area A is raised to 
potential @ and is placed at a distance d from a similar parallel plate 
which is earthed. Neglecting edge effects, show that it is attracted to the 
earthed plate with a force e, Ad?/2d. 


Solution. Taking the origin O on the earthed plate and the x-axis 
perpendicular to this plate (Fig. 1.11), if we neglect distortion in the region 
of the edges, the potential Vwill depend upon x alone, Thus, between the 
plates, V will satisfy V?V = d?V/dx* = 0 and, hence, V = A+Bx. Since 


x 


Veo 
Fig. 1.11 
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V =Oover x = 0, V = ¢ over x = d, these boundary conditions require 
that A = 0, B = ¢/d. It follows that 

V = x/d. (1.38) 
This result shows that E is parallel to Ox and that its magnitude is 
E = —dV/dx = —#/d. It follows from equation (1.21) that the surface 
density of the charge on the side of the insulated plate facing the earthed 
plate is given by ¢ = ¢, ¢/d and this corresponds to a force per unit area 
of o?/2&) = && #?/2d. Thus, the total force applied to the insulated plate 
is 6) Ad?/2d. o 


Problem1.15 A conducting sphere is cut into two hemispheres, earthed, 
and placed in a uniform field with the cut normal to the field, Calculate 
the force tending to separate the hemispheres. 


Solution. The field has been calculated in Problem 1.12. The surface 
density of charge induced on the sphere was shown to be 3é Ey cos 0. 
Thus the force per unit area applied to the surface is 34) EZcos”0. 

Divide one of the hemispherical surfaces into ring-shaped elements by 
planes perpendicular to its axis of symmetry as indicated in Fig. 1.12. 

$ cos 


AN 


Then ad is the width of a ring and 2za sin 0 is its circumference. Hence, 
its area is 2na?sin 0 d0. The resultant force on the ring must be along the 
axis of symmetry @ = 0 and therefore only the component of the pressure 
eq E2cos?0 in this direction, namely $2 Ejcos*0, contributes. Thus, the 
resultant force on the ring is 9¢) a*E3cos*0 sin 0 d0. Integrating for all 
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rings belonging to the hemisphere, we obtain for the force tending to 
separate the hemispheres 


[2 90 na*EZcos°0 sin dO = ay EBna?. fal 


Problem 1.16 The centre of the inner sphere of the system described in 
Problem 1.11 is displaced a small distance 5. Calculate the resultant force 
on this sphere to the first order in 6. 


Solution. Let O, O' be the centres of the inner and outer spheres respec- 
tively. Let (r,0) be spherical polar coordinates relative to O as pole and 
OO’ as reference line. In the region between the spheres, we assume that 
the potential is given by 


v= A+2+( crt cov. (1.39) 


When 6 = 0, this must reduce to the form of equation (1.26); C and D 
must accordingly be of order 5. Over r = a, V = ¢ for all 0; this boundary 
condition leads to the equations A+B/a = ¢, Ca+Dj/a? = 0. On the 
outer sphere, r = b+écos0 to the first order in 6 and V must vanish. 
Substituting in equation (1.39) and approximating to the first order in 5, 
we derive the condition 
B 6 D 
A+F (1-§ cose) +(co-+j2 ond =0. 
This is satisfied for all 0 provided A+B/b = 0, Ch+D/b? = 5B/b*. A, B, 
C, D can now be found. A and B have the values already calculated in 
Problem 1.11 and aC = —D = (da*bd)/[(b—a)(b? —a°)]. 
The charge density on the inner sphere is next calculated to be given by 


r=) -oe(2- 


The resultant force on the sphere is now determined as in the previous 
problem to be 


) os cos 0. 2na?sin 0 dO = 2me (= -c) f cos"0 sin 0 dO 
0 28 ¢ ) 
_ 4 Mtq 6a*b? op? o 
© (bay (b=) 


1.6 Images In some problems it is possible to imagine that the field is 
generated by point charges alone and to show that this field satisfies all 
the boundary conditions. Some of the point charges used will have no 


19 


real existence and are called images; nevertheless, their combined potential 
will satisfy Laplace’s equation and, provided it also satisfies the necessary 
boundary conditions, it follows from the uniqueness theorem that it 
represents the correct solution to the problem. 

A point charge q placed at a point A distant d from an infinite earthed 
conducting plane provides a simple example. The potential generated 
by this charge and an image charge —q situated at the optical image A’ 
of A in the plane (i) vanishes on the plane, (ii) vanishes at infinity, (iii) posses- 
ses the appropriate singularity in the neighbourhood of the charge at A. 
By the uniqueness theorem, it is accordingly the correct potential for the 
field between the charge and the plane. 


Problem 1.17 Two semi-infinite, earthed, conducting planes meet at 
right angles at a common edge. A charge q is placed at a point A ona 
bisector of the angle between the planes at a distance d from the common 
edge. Show that the charge is attracted towards this edge with a force 
(2/2—1)q?/(L6ne d?). 


Solution. The image system for the field in the angle between the planes 
comprises charges —4q, q, —q at points B,C, D respectively, where ABCD 
is a square (Fig. 1.13). This system of charges clearly generates zero 
potential over both planes as required. Thus, if the planes were removed 
and the image charges placed in position, the field in the neighbourhood 


Fig. 1.13 
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of A would be unaltered. It follows that the force on the charge at A would 
be unchanged. This force can now be easily calculated by application of 
Coulomb’s law and the result stated derived. 


Problem 1.18 A point charge q is placed outside an insulated, uncharged 
sphere of radius a at a distance f from the centre. Calculate the surface 
density of the charge on the sphere at a point distant r from q. 


Solution. The image system comprises a charge —aq/f at the point A’ 
inverse to A in the sphere and a charge aq/ fat the centre O (Fig. 1.14). For 
the potential generated by the three charges at a point P on the sphere is 


1 (gaa 1 ag 

a(t Prt ya} (1.40) 
Now, OA.OA’ = OP?; hence, the triangles OAP, OPA’ have the angle 
at O in common and the sides about this angle in the same ratio, ie, 


Fig. 1.14 


OA/OP = OP/OA’. These triangles are similar, therefore, and this 
implies that r/r’ = f/a. The expression (1.40) accordingly reduces to 
q/(4ne, f), showing that the three charges yield a constant potential over 
the sphere. Further, since the image charges enclosed by the sphere cancel, 
it follows from Gauss’s law that the flux of electric intensity out of the 
sphere generated by the three charges is zero. But, if is the surface density 
of charge on the sphere S, then 


charge on sphere = j ao dS = & { E, dS = 0. 
Ss Ss 
Thus, all the boundary conditions are satisfied by the image system. 
Resolving the intensities at P due to the charges along OP, we find that 
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1 faq 1 aq 1 q 
E, = eft SF reoeo—Somo] (1.41) 


Substituting r’ = ar/f, cos@ = (r? + f?—a*)/2rf, 
cos $ = (f?—a—r*)/2ar and noting that ¢ = &>E,, equation (1.41) 


gives 
iy ES 
-£[5-S | o 


Problem 1.19 A dipole of moment m is placed at a distance f from the 
centre of an earthed sphere of radius a. The axis of the dipole points 
directly away from the sphere. Show that the charge induced is am/f “ 


Solution. Suppose the charges —q, q of the dipole lie at the points A, 
B respectively (Fig. 1.15), where OA = f, OB =f+é6f. Then m= qof 
and 5f — 0. The inverse points A’, B’ are such that 
a , @ =( *f) 
5 B’ = ——_—._ = —( 1-—— 
ii, A Geey. Saad 
to the first order in 5f. Hence, A’B' = a5 f/f?. The sphere will be at 
zero potential if we introduce image charges aq/f at A’ and —aq/(f +f) 
at B’. To the first order in df, —aq/(f+6f) = —(aq/ f)+(aq/f?)df. Thus, 
the charges at A’, B’ are equivalent to a dipole of moment 


o 


OA’ = 


Fig. 1.15 


aq.A'B'/f = a'q5f/f? = a°m/f® directed along B'A’ and a further 
charge aq5f/f? = am/f? at B’. In the limit, as 5f — 0, both the image 
dipole a°m/f* and the image charge am/f? are situated at A’. 
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The charge induced on the sphere is the total image charge enclosed, 
ie. am/f?. Oo 
Problem 1.20 A closed hemispherical shell of radius a is constructed 
from sheet metal and is earthed. A charge q is placed inside the shell, on 
its axis, a distance f from its centre. Find the force experienced by the 
charge. 


Solution, If q is placed at A, image charges —aq/f, —q, aq/f are intro- 
duced at the points A’, B, B’ respectively (Fig. 1.16). A, A’ are inverse 


Ky -aq/t 


Blagft Fig. 1.16 

points with respect to the completed sphere; B is the optical image of A 
in the plane face of the hemisphere; B’ is the inverse of B with respect to 
the sphere. Clearly, B’ is the optical image of A’ in the plane face. The 
charges at A and A’ and the charges at B and B’ make the potential on the 
curved surface of the hemisphere vanish. Also, the charges at A and B and 
the charges at A’ and B’ make the potential zero over the plane face, Thus, 
the image system solves the problem. . 

In particular, the resultant force on the original charge at A due to the 


images is 
ef Ly er rafiota) valh | 
Snel 4f? (frase F-a/fF 
directed towards B. 
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Problem 1.21 An infinite line charge of density q per unit length is fixed 
parallel to the axis of an infinite circular conducting cylinder of radius a. 
The distance between the line charge and cylinder axis is f. Calculate the 
charge per unit length of the cylinder if no force is exerted on the line 
charge. 


Solution. A right-section of the cylinder and line charge is shown in 
Fig. 1.17. The line charge and cylinder axis meet the plane of section in A 


and O respectively. Introduce parallel image line charges —q, Q+q to 
meet the plane of section in A’ and O respectively, A’ being the point 
inverse to A in the circular section. Then, at any point P on the cylinder, 
the potential generated is given by 


= —2q log AP +2q log A’P—2(Q+q)log OP 
AP 

= 2aloe( ip) 10+ aoga 

= —2qlog f —2Q loga. 


This is constant over the whole surface of the cylinder as required. Since 
the total image charge per unit length enclosed by the cylinder is Q, this 
is its charge per unit length. 
The field intensity at A due to the two images is 

ek [2072 ca } 

4ney fF - eS 
This vanishes if Q = a°q/(f?—a”). oO 
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Problem 1.22 A pair of infinite circular conducting cylinders, both of 
radius a, are placed with their axes parallel a distance 2d apart. Their 
potentials are +@ and —@. Calculate their charges. 

Solution. A right-section of the cylinders by a plane meets their axes in 
Oand O ' (Fig. 1.18). A, A’ are points on OO’ which are inverse with respect 
to both circular sections. Thus, if OA’ = O'A = x, then 


P 


(a hee 


Fig. 1.18 


OA = O‘A' = 2d—x and, hence, x(2d—x) =a’. Solving this last 
equation for x, we find x = d+,/(d?—a’). (N.B. positive root since 
x>d) 

Introduce parallel infinite image line charges q, —q to intersect the 
plane of section at A, A’ respectively. The potential generated at any point 
P on the left-hand cylinder is now given by 


Pie | q oe | x 
= ag, WERE te Ree sie (*). 
This is constant, as required. To give the correct potential, it is necessary 
that 
= 2ne 
4 iogid/a+ /(@?/a*—1))" 

The potential over the other cylinder can now be shown to be —#, 
similarly. 

The last equation gives the charge per unit length over the left-hand 
cylinder. The charge on the other cylinder is —q per unit length. 


1.7 Capacitors. Energy An insulated conductor specially designed to 
store electric charge is called a condenser or a capacitor. Such a device 
usually comprises two plates, one being earthed and the other being 
insulated and serving to hold the charge. It is found that the charge Q 
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stored is directly proportional to the potential @ above earth of the 
insulated conductor; the ratio Q/¢ is called the capacitance. The purpose 
of the earthed conductor is to increase the capacitance. 

Making use of results stated in § 1.2, Problems 1.11 and 1.14, the capaci- 
tances of (a)a conducting sphere of radius a, (b) the same sphere surrounded 
by an earthed spherical shell of radius b, and (c) a parallel plate capacitor 
whose plates have area A and are separated by a distance d, are all easily 
found. The results are (a) 47&9 a, (b) 472 ab/(b—a), (c) & A/d. 


Problem1.23 The axis of an infinite circular conducting cylinder of 
radius a is parallel to an infinite earthed plane and distant d from it. 
Find the capacitance per unit length of the cylinder. 


Solution. The image line charges g, —q at A and A’ in Fig. 1.18 have 
been shown (previous problem) to generate a constant potential @ over 
the left-hand cylinder, provided 

2nig 
= Toold/at(d2/ae—1)F1~ 
41 jos{d/a+(d?/a?—1))] 
They also generate zero potential over the plane which is the perpendicular 
bisector of OO’. This image system accordingly solves our present problem. 

By Gauss’s law, q will be the charge per unit length on the cylinder. 

Hence, its capacitance per unit length is 


ee 2n& 
@ ~ log(da+(@/a—1)" 5 


If charges q,,42,...,4, are situated at points P,,P,,...,P, and 
is the potential at P; due to all charges except q,, the potential energy U of 
the system is given by 


u=4¥ aM. (1.42) 


This is the work which would be done by the attractions and repulsions 
between the charges if these were separated to be at great distances from 
one another. 

In the case of a capacitor C whose insulated plate has been raised to 
potential #, the charge on this plate is Cd; the contribution of this charge 
to the energy of the system is therefore }C#*. The earthed plate is at zero 
potential and the induced charge upon it accordingly makes zero contri- 
bution. Thus, the total energy of the capacitor is $C?. 


Problem 1.24 Calculate the internal electrical energy of a proton 
assuming that its charge ¢ is distributed uniformly over a sphere of radius a, 
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Solution. Employing the result found in the solution to Problem 1.9, the 
field intensity in the proton will be radial and of magnitude er/(4ze9 a’). 
Since E = —dV/dr, this implies that V = C—er?/(8n&9 a*), where C is 
constant. Outside the proton, the field is the same as if all the charge were 
at the centre; thus, V = ¢/(4neqr). Since V must be continuous over the 
surface r = a, C = 3e/(87& a). 

Consider the charge enclosed between spheres of radii r, r+dr with 
their centres at the centre of the proton. The quantity of charge is 3er2dr/a* 
and its potential is e(3a?—r?)/(8n¢, a*). Its contribution to the internal 
energy is therefore 3¢?r?(3a—r?)dr/(16meqa°). Integrating over the 
interval 0 < r < a, the total internal energy is found to be 3e*/(208 a). 


Problem 1.25 A spherical conductor of radius a is kept by batteries at a 
constant potential ¢, It is surrounded by a concentric, spherical conductor 
which is earthed. Under the internal forces this latter conductor contracts 
from radius b to by. Calculate the work done by these forces. 


Solution. The charge on the inner sphere has already been calculated 
at equation (1.29). Thus, the energy of the system before contraction is 
28 abp?/(b—a). 

After contraction, the energy is 2&9 ab, ¢?/(b,—a) and the gain in 
energy is accordingly U = [2neo a2(b—b,)o*|/[(b—a)(b, —a)]. 

‘The increase in charge on the inner sphere is 


Arey abs > _4neyabp _ 4m a(b—bi\p 
b,-a b-a (b—aXb,—a@) * 


This charge is raised by the batteries from earth potential to potential #. 
Thus, the work done by the batteries on the system is 2U. Of this, one half 
appears as a gain in energy of the system and the other half corresponds to 
the work done by the internal forces. Hence, the work done by the internal 
forces is U as given above. Oo 


EXERCISES 


1. A line doublet comprises a pair of equal and opposite parallel infinite 
line charges a small distance d apart. If q and —q are the charges per unit 
length, show, by neglecting powers of d after the first, that the potential 
of the field is given by V= mcos $/(2nzo p), where m = qd is the moment 
of the line doublet and (p, ¢) are cylindrical polar coordinates in a frame 
Oxyz for which Oz is parallel to the line charges and midway between 
them and Ox lies in the plane of these charges. 
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2. Prove that the field intensity due to a uniformly charged ring of total 
charge q and radius a, at a point on its axis distant z from its centre, is 
qz(a? +27)" 3/(4neo). 

A fixed cylindrical shell of radius a and length 2I carries a total charge 
Q uniformly distributed over it. A particle of mass m and charge —q is 
free to perform small oscillations about the centre of the cylinder’s axis, 
along the line of this axis. Prove that the period is 2x[4e m(a? + [7)3/Qq]. 


3. Infinite coaxial conducting cylinders have radii a and b (a < b). The 
outer cylinder is earthed and the inner is raised to a potential ¥,. Assuming 
that the potential between the cylinders takes the form shown in equation 
(1.22), determine the field. Hence show that the capacitance per unit 
length of this cylindrical condenser is 27#9/[log(b/a)]. 


4. If the conducting sphere in Problem 1.12 is insulated and given a 
charge Q before being placed in the uniform field Ey, show that the poten- 
tial of the resultant field outside the sphere is given by 


Ve e(3- t)poon0e 2 


4negr 
Deduce that the surface density of the charge on the sphere is positive 
everywhere provided Q > 12xa7% Ey. (Assume Ey > 0). 


5. A conductor in the form of an infinite plane with a hemispherical 
boss of radius a, is earthed in the presence of a point charge q at a point P. 
The line OP makes an angle @ with the normal to the plane, where O is 
the centre of the circular base of the hemisphere and OP = f- Show that, 
if a/f is small, the component of attraction on the charge normal to the 
plane is greater than it would be in the absence of the boss by an amount 
of approximately q?a°cos 6(1 +3 cos 26)/f*. 


6. The centre of an insulated, uncharged conducting sphere of radius 
a is at the midpoint of the straight line joining two equal point charges, 
which are a distance 2f apart. If a/f is small, show that the force on either 
charge would be increased in the approximate ratio (14+24a5/f):1, if 
the sphere were removed. 


Chapter 2 
Polarized Media 


2.1 Poisson’s Distribution When an insulator or dielectric is placed in 
an electric field it becomes polarized, i.e. any small element behaves like 
a dipole and contributes to the resultant field. This effect is due either to a 
slight separation of the positive and negative charges within the atoms or 
to the alignment of the axes of existing atomic dipoles under the influence 
of the field. 

If dv is any small element of dielectric volume and P dv is its dipole 
moment, then P is called the polarization of the dielectric at the position 
of the element. 

Ifa dielectric is bounded by a surface S, the field at any external point 
due to the polarization within S is identical with that due to a surface 
charge over S having density Pn (n is the unit outwardly directed normal 
to S) and a charge distributed over the interior of S with volume density 
—div P. This is Poisson's distribution. 

Inside a dielectric, E is defined to be the field intensity of the Poisson 
distribution. The actual force exerted upon unit charge situated inside 
a dielectric is difficult to define unambiguously. If the charge is supposed 
to be placed in a small cavity, the force depends on the shape of the cavity. 
If the cavity is needle-shaped, with its axis parallel to P, the force can be 
shown to be E as just defined. 


Problem 2.1 Calculate the field due to a uniformly polarized sphere. 


Solution, Taking the centre of the sphere O as the origin for spherical 
polar coordinates (r,0), we shall suppose that the polarization P is 
always directed along the z-axis (Fig. 2.1). 

Then, at any point (r,0,) on the sphere’s surface, P.n = Pcos0; 
this gives Poisson's surface distribution. div P = 0, so that there is no 
volume distribution. 

The dielectric sphere can now be imagined to be removed and replaced 
by the surface distribution; the methods of Chapter 1 are then applicable. 
Let V, be the field potential within the sphere and V, the potential outside. 
Both potentials must satisfy Laplace’s equation and we shall take 


B 
V, = Arcos 6, Vz = 720080. 


Fig. 2.1 


The boundary conditions to be satisfied over r = a are (i) V is continuous 
hence Aa = B/a? and (ii) E,2—E,; = o/% = P.cos 0/éo ; since 
E, = —@V/ér, the second condition yields A+2B/a* = P/e. Thus, 
A = P/3eq and B = Pa®/36. 

These results imply that the internal field is uniform and of intensity 
—P/3eq and that the external field is identical with that due to a dipole 
of moment 4xPa*/3 placed at O. 


2.2 Fields in Dielectrics In most dielectrics, it is found that P is directly 
proportional to E. Thus, we write 

P = xeoE, (2.1) 
where is called the electric susceptibility. 

Replacing all polarized dielectrics by their Poisson distributions, the 
situation reverts to that studied in Chapter 1 and the field equations used 
in the absence of dielectrics become applicable. In particular, 

E = —grad V. (22) 
Also, at any point there will be a density p of real charge and a density 
—div P of imaginary Poisson charge. It follows from equation (1.17) that 
div E = (p—div P)/e,. This equation can now be written in the form 
divD = p, (2.3) 
where D = 6 E+P = 6(1+7)E, by equation (2.1). Writing e = 1+y, 
this gives 
D = cE. (2.4) 
& is called the dielectric constant. 
Equations (2.2)-(2.4) are the fundamental equations of the electrostatic 
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field in the presence of dielectrics. D is called the electric displacement. If a 
small disc-shaped cavity is cut in a dielectric with its axis aligned with P, 
the field intensity inside is D/ég . 

Problem 2.2 The potential for a field in a uniform dielectric whose 
constant is ¢, is given by V = V(r), where r is distance measured from a 
point O. Calculate the charge density p. If p = po (a/r), find V. 


Solution. Equation (2.2) gives E = —V'r/r, where V’ = dV/dr. By 
equation (2.4), D = — 2&9 V'r/r and it then follows from equation (2.3) that 


Pe ss withpiia: ew enoun 
p= 1% ao“) = —at9(V"+2V'r) = —F (PV), 
In the case p = po(a/r)’, integrating the last equation we find that 


where A is constant. It is easy to verify that the charge inside a small sphere 
of radius 5 and centre O tends to zero with 6, i.e. there is no point charge 
at O. Hence A = 0. A further integration now gives 


Vv VED tc citait, Oo 
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Equation (2.3) implies that the flux of D out of any closed surface $ 
is equal to the total (actual) charge enclosed in S. This is the generalized 
form of Gauss’s flux law (equation (1.16). 


Problem 2.3 Amend Coulomb's law (equation (1.3)) for the case when 
the two charges are immersed in an infinite dielectric of constant ¢. 


Solution. Consider the situation when the charge q is alone in the dielec- 
tric at a point O. The resulting field will be spherically symmetric about O 
and directed radially from this point. Thus, the flux of D out of a sphere 
of radius r and centre O is 4nr?D. By Gauss’s law, 4ar7D = g. Hence, at a 
point P where OP = r, 


= ior (23) 
Employing equation (2.4), this gives 
ae | 
E= ann” (2.6) 
Equation (2.2) then leads to 
ea | 
V= aaa (2.7) 
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Now introduce the charge q’ at P (assumed placed in a radial needle- 
shaped cavity so that it is free to move in the direction of the applied force). 
The force exerted upon it is given by 


44 
4nee r° 
This is the amended Coulomb's Law. o 


F=qE= " (2.8) 


2.3 Harmonic Functions Eliminating E and D between equations 
(2.2)+(2.4), we find that 


div(e grad V) = —p/eo. (2.9) 
Over any region in which « does not vary, this reduces to 
VV = —plety. (2.10) 


Thus, if p vanishes over such a region, V will be a harmonic function. 

Ifa surface S carries a charge having density o and the regions on either 
side of S are occupied by dielectrics of constants 2, ,£2, the following 
boundary condition is applicable: 

Dy,—D,, = 4. (2.11) 
This follows from Gauss’s law by an argument similar to that used to 
derive equation (1.20), The subscript n indicates a component along the 
normal to S directed from the region 1 to the region 2; this is the amended 
form of the condition (1.20). If S is the surface of a conductor, then D, = 0 
inside the conductor and D, is directed along the normal; thus, the condi- 
tion reduces to 
D,=<a. (2.12) 
The potential V is assumed to be continuous across S; this is equivalent 
to the statement that all tangenitial components of E are continuous. 


Problem 2.4 Two conducting closed surfaces, one inside the other, form 
the plates of a capacitor. If the region between the surfaces is filled with 
uniform dielectric of constant ¢, show that the capacitance is multiplied by 
8 


Solution. Let Sy be the outer surface and S, the inner surface. With 
no dielectric between the plates, suppose So is earthed and S, is raised to 
potential . Then the potential V in the region between S, and S, is 
harmonic and satisfies the boundary conditions V = 0 on Sy and V = 
on S,. These conditions serve to determine V uniquely. When V has been 
found, the surface density of the charge on S, can be found from the 
equation ¢ = & E, = —&0V/dn, 2/dn denoting differentiation along the 
outward normal to S,. 
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Now suppose the dielectric is introduced and the plate S, is maintained 
at the potential }. Then V satisfies the same conditions as before and is 
therefore unchanged. However, the surface density of the charge on S, 
is now given by o’ = D, = —et9 0V/On = ea, ie. the total charge on S, 
is multiplied by «. Since the potential difference between the plates is still 
¢, this means that the capacitance has been increased by a factor e, o 


Problem 2.5 A conducting sphere of radius a is coated with a dielectric 
shell of thickness d and constant «. It is then surrounded by a concentric 
conducting shell of radius b, which is earthed. Find the capacitance of the 
inner sphere. 


Solution. Suppose the inner conductor is raised to potential $. In region 
1 (Fig. 2.2), we assume V = A+B/r and in region 2 we take V = P+Q/r. 
Over r = a, V = and over r = b, V = 0. Also, V is continuous across 
r=(a+d). These conditions give ¢ = A+(B/a), 0 = P+(Q/b), 
A+[BMa+d)] = P+[O/a+d)]. 


Fig, 2.2 =< 


In region 1, the field is radial and, hence, E = —dV/dr = B/r?; thus, 
D = e&B/r’. Similarly, in region 2, D = &9Q/r?. Applying condition 
(2.11) to the surface of discontinuity r = a+d, we get 
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We can now solve for A, B, P, Q thus: 


n peaea= ay, yiiee ab(a+d) ¢ 
bd+ae(b—a—d)"* bd+asb—a—d)"* 
P = —eB/b, Q = cB. 

Using equation (2.12), we calculate that the density of charge on the 
sphere r =a is eg) B/a* and, hence, that the total charge is 47ee, B. 
Thus, the capacitance is 

AnetgB Aree abla+d) 
@ — bd+aeb—a—d)” 

A special case is when the dielectric fills the gap between the plates. 
Then, d = b—a and the capacitance is 4mee, ab/(b—a). Comparison of this 
result with that given in §1.7 reveals that the presence of the dielectric 
again increases the capacitance by a factor e. o 


A 


Problem 2.6 A dielectric sphere of radius a and constant ¢ is placed in a 
uniform field of strength E,. Find the resultant field produced. 


Fig. 2.3 


Solution. Taking the origin of spherical polar coordinates (r,0,) at 
the centre of the sphere and z-axis in the direction of the applied field 
(Fig. 2.3), in the absence of the sphere the potential is given by 
V = —E,z = —Eyrcos@. When the sphere is introduced, its polariza- 
tion will distort the uniform field but the distortion will be negligible at 
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great distances; accordingly we shall assume that, outside the sphere 
(region 2), 


V=h= —Eyrcos0+4 cos0. 


This function is axially symmetric about Oz and satisfies Laplace's 
equation throughout region 2. 

Inside the sphere (region 1), the potential must be axially symmetric 
and finite at r = 0. Hence, we assume V = V, = Brcos0. 

Across the surface r = a, both V and D, must be continuous. Since 
D, = tE, = —s&o(dV/dr),-,, these conditions are satisfied provided 

—E,a+(A/a?) = Ba, Ey+(2A/a*) = —eB. Thus, A = a°Eo(2—1)/(e+2) 
and B = —3Eo/(e+2). These give 
e—1 aE, 
5 ie —Eorcos+ 53 008 0. 

The first of these results implies that the field inside the sphere is uniform 
of intensity E = 3E,/(e+2). The second result implies that the polarized 
sphere contributes to the external field like a dipole at O of moment 
Ane, a®E,(e—1)(e+2). im} 


Problem 2.7 A dipole of moment m is placed at the centre of a spherical 
cavity in an infinite dielectric of constant e. Show that the field intensity 
in the dielectric is the same as that produced by a dipole of moment 
3m/(2e+ 1) in empty space. 

Solution. Take spherical polar coordinates (r, 6) as indicated in Fig. 2.4. 
Inside the cavity, the potential will be the sum of a part due to the dipole 
and a part due to the polarized dielectric. We assume 
(mcos 0) 
(rt r°) 
In the dielectric, the potential must vanish at infinity and we shall assume 
V = Vz = Bcos6/r*. 

Across the wall of the cavity r = a, V and D, are continuous. These 
conditions lead to the equations 

m 


m B 
wae +Aa= 3, 5 aa bear 5 
Hence, A = —m(e—1)/[4ze9 a*(2e+1)], B = 3m/[4z60(2e+1)]. The po- 
tential in the dielectric is therefore given by 
3m cos 0 
* GregQe+ DF 


V=,= + Arcos 0. 


yy 
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Fig. 2.4 

and this is the potential of a dipole of moment 3m/(2e+ 1) placed at O in 
empty space, Oo 
2.4 Images 


Problem 2.8 The semi-infinite region on one side of a plane is occupied 
by dielectric of constant «. A point charge q is placed outside the dielectric 
at a distance d from the plane. Find the force on q. 


Fig. 2.5 
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Solution. If q is at A (Fig. 2.5), let A’ be the optical image of A in the plane 
face of the dielectric, We shall assume that the field potential outside the 
dielectric can be generated by point charges q at A and q’ at A’. Both 
charges are to be thought of as situated in empty space. Thus, if r, r’ are 
distances measured from A, A’ respectively, then 


V= am (+5) r>r (2.13) 


Inside the dielectric, V will be supposed generated by an image charge 
q; at A (also regarded as being alone in space). Thus V = q,/(47&9 r), 1’ < r. 
V must be continuous at the plane face, and since r = r at Q, 
q+d =4%- (2.14) 
Outside the dielectric at Q, the component of E in the direction of the 
normal n due to the charges gq, q' is 
q q 
Rare deat cos 0. 
Multiplying by &9, this gives D,. Inside the dielectric at Q, 
E, = —q,00s 0/4ne)r?, Hence, D, = —2q, cos 0/4nr*. But D, is con- 
tinuous at Q and, therefore, 
qq = —8y. (2.15) 
Solving equations (2.14), (2.15), we find that 


2 & 
Gh se Ce (2.16) 


By appeal to a uniqueness theorem, it now follows that the fields in the 
two regions are correctly generated by these two systems of image charges. 

In particular, if the dielectric is removed and the image q/ is introduced, 
the field outside the dielectric will be unaltered. Hence, the force on q is 
equal to the attraction due to q’, namely (e—1)q?/[16neo(e+1)d?]. O 


Problem 2.9 If, in the previous problem, a real point charge e is placed 
in the dielectric at the point A’, calculate the force on q. 


Solution. We suppose the potential outside the dielectric to be generated 
by point charges q, q’ at A, A’ respectively, as before. Inside the dielectric 
we suppose the potential generated by charges q,, e/e at A, A’ respectively; 
it is understood that all charges are being imagined situated in a vacuum, 
so that the potential arising from this pair of charges is 


ch Cheng 
4negr 4negyr’? (2.17) 


av 


comparing the second term with equation (2.7), it will be seen that the 
behaviour of V in the neighbourhood of the charge e at A’ has been 
generated correctly. 

We now require that the potential functions (2.13), (2.17) should be 
continuous across the plane r = r’ and, further, that the associated D, 
should be continuous, These conditions lead to the equations 

qt+q =qte/s, 4'-q =e—a,. 

Solving for q’, q, , we find that 
a wee pte us 
rs heres Ss Rr ba 
The force on q is the attraction of q’, ice. 

aq _ a(e~1)q~2] o 

16nd? 16me) (e+) 

Problem 2.10 An infinite line charge q is placed parallel to the axis of an 
infinite circular cylindrical dielectric of radius a. If the distance between the 
line charge and the cylinder axis is f(> a), calculate the force per unit 
length on the charge. 


{b) 


Fig.26 
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Solution. It will be assumed that the potential outside the cylinder can 
be generated by the actual line charge g through the point A (Fig. 2.6(a)) 
and parallel image line charges (i) g’ through the point A’ inverse to A 
in the cylinder and (ii) —q’ through O along the cylinder axis, For the 
actual field, the flux of D out of any length of the cylinder must be zero, 
since it contains no charge. The overall image charge within the cylinder 
must be taken to be zero, therefore, as we have assumed. The potential 
inside the cylinder will be supposed to be generated by a line charge q, 
through A (Fig. 2.6(b)). As usual, this charge is thought of as being alone 
in space. 

The potential must be continuous at all points P on the cylinder. At such 
a point just outside the cylinder, the potential is being assumed to be 
(—qlogr—q'logr’+q' loga+C)/2ne,. The potential due to a line charge 
is always arbitrary to the extent of an added constant and advantage has 
been taken of this to introduce the constant C. The potential inside the 
cylinder at P is being assumed to be —q, logr/2xe. The condition for 
continuity is accordingly (q—q,)logr+gq'logr’ = q'loga+C. But, by 
similar triangles, r’ = ar/f and, thus, this condition can be written 
(q—4:+q'logr = q'logf+C. If this is to be satisfied for all r, it is neces- 
sary to take 


u-F =4q (2.18) 
and C = —q'logf. 
Continuity of D, at P requires that 
4 Le ae ety 
Faq 080-7 008 H+5— = 5080. (2.19) 


By projecting OA on to AP, it follows from the figure that 
fcos@ = acos6+r. Since f = ar/r’, this implies that 
1 1 1 
7oOsd = —cosO+=. 
Equation (2.19) according reduces to 


F cos 0 
(q—q'—e4,) 5 = 0. 
This is satisfied at all points on the cylinder, provided 


+e, = 4. (2.20) 
Solving equations (2.18), (2.20), we find that 


39 


The field intensity at A due to the image charges within the cylinder is 
now found to be 
gel @ 
2nty e+1 f(f?—a") 
directed towards O, Thus, the force per unit length of the line charge at 
Ais a7q*(e—1)/[2neo(e+ 1) f(f? —a?)]. o 


2.5 Magnetostatic Field The sources of the magnetostatic field are 
certain atoms which behave towards one another in the same manner as 
electric dipoles. However, the physical basis of this magnetic property is 
quite different from that of an electric dipole; it is caused by the circulation 
of electric charge within the atom. Nevertheless, the similarity of behaviour 
is so close that the mathematical theory we have developed to describe the 
phenomenon of electric polarization can be applied with only small 
modification to magnetically polarized materials. 

Although isolated magnetic poles have never been observed, it is 
mathematically convenient to base the theory on a law of force between 
two poles of strengths p,, p2. This is the analogue of equation (1.3) and is 
written 

Ho Pi P2 
F= 4a? f, (2.21) 
where jig = 1:257 x 10~° is called the permeability of free space and p,, p2 
are measured in SI units (see p. 95).* 

In free space, the force per unit pole is denoted by B and is termed the 
magnetic induction, Thus, for a hypothetical pole of strength p at the 
origin, it follows from equation (2.21) that B = jg p?/4nr?. Also, the 
couple which acts upon a magnetic dipole of moment m placed in a field 
BismxB. 

In free space, the magnetic intensity H is defined to be B/jig. Thus, for a 
pole of strength p at the origin, H = p?/4nr?. As has been remarked in 
§1.2, this particular field is conservative and can be derived from the mag- 
netic potential 2 = p/4zr by application of the equation 


H = —gradQ. (2.22) 


Since any magnetostatic field can be generated by a distribution of poles 
(in dipole pairs), equation (2.22) is certainly valid at all points in free space. 
In particular, for an isolated dipole of moment m at the origin, 


* ‘9 is placed in the numerator to permit B to be taken as the force per unit pole instead 
of H. This is the accepted modern convention. 
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Comparing this equation with equation (1.9), it should be observed that 
the potential of an electrostatic distribution can always be converted into 
the potential of the corresponding distribution of magnetic poles by 
setting 9 = 1. 

It can now be shown, exactly as in the corresponding electrostatic 
situation that, if M is the magnetic moment per unit volume within a 
magnetic material, the field generated externally is identical with that of 
a distribution of poles of volume density — div M over the region occupied 
by the material and a distribution of poles of density M.n over its surface 
(n is the unit outwards normal). This is the Poisson distribution for the 
material. M is termed the magnetization. Within the material, the potential 
Q and intensity H are defined to be those of the Poisson distribution. This 
ensures that equation (2.22) is valid everywhere. 

Since H = p?/4nr? for a pole p at the origin, the flux of H out of any 
closed surface in a field due to a hypothetical distribution of poles is equal 
to the total magnetic charge within the surface. As shown in §1.4, the 
differential form of this result is divH = p, where p is the density of 
magnetic charge. Application of this result to the Poisson distribution 
leads to the equation divH = —div M; hence, 


divB =0 (2.24) 
where B= u(H+™). (2,25) 


This definition of B within a magnetic material is consistent with the 
equation B = oH in free space. 

It can be proved that H is the intensity within a small needle-shaped 
cavity cut inside the material with its axis parallel to M. Thus, the force 
on unit pole in such a cavity is 49 H (not B). The intensity inside a small 
penny-shaped cavity whose axis is aligned with M is found to be H+M. 
Thus, the force on unit pole in the cavity is zo(H+M) = B. 

Equation (2.24) is the differential form of Gauss’s law for magneto- 
static fields; this states that the flux of B out of any closed surface (which 
may intersect magnetized materials) always vanishes. 


Problem2.11 A sphere with centre at the origin and radius a is magnetized 
so that the components of magnetization are A(x, y,z+a). A dipole of 
moment m is placed at the point (c,0,0) (c > a) with its axis parallel to 
Ox, Calculate the couple acting upon it. 
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Solution. Since M = A(x,y,z+a), divM = 3A; thus the Poisson 
distribution over the interior of the sphere has density — 3A. If r = (x,y,z) 
is the position vector of a point, M = A(r+ak), where k is the unit vector 
along Oz. Thus, the radial component of M is M.r/r = A(r+acos8), 
where @is the spherical polar angle between r and k. Hence, on the sphere’s 
surface, M,, = 2a(1+cos@) and this defines the surface Poisson distri- 
bution. 

The volume charge of density — 3A and the part of the surface charge of 
density 4a can, for the purpose of calculating the external field, be con- 
densed to a point charge at O of magnitude $na*(—34)+4na*(Aa) = 0. 
It follows that the external field is that generated by a surface charge of 
density Aa cos 0. 

The potential of a spherical surface distribution of this type has already 
been found in Problem 2.1; setting P = Ja and gq = 1 in the result of 
this calculation, we get 2 = Aa*cos 6/3r? as the external potential. Thus, 
the external field is the same as that due to a dipole of moment 47/a*/3 at 
O with its axis along Oz. 

The radial and transverse components of H due to this dipole can now 
be found by calculating —62/ér and —02/r 00. In particular, the intensity 
at (c, 0, 0) is found to be Aa*/3c? directed along zO. Thus, B = —kApy a*/3c* 
at this point and the given dipole experiences a couple 
mxB = mix B = jmdyoa*/3c (i and j are unit vectors along the 
x- and y-axes. oa 


2.6 Harmonic Functions For most materials, M is directly proportional 
to H and we can write 
M = 4H, (2.26) 
where x is the magnetic susceptibility (c.f. equation (2.1)). Using equation 
(2.25), this leads to the result 
B= po(l+2H = pyoH, (2.27) 
where = 1+ 7 is termed the permeability of the material. 
In such materials, the fundamental equations of the magnetic field are 
H=-—gradQ, divB=0, B= pyoH. (2.28) 
Q therefore satisfies the equation 
div(u grad Q) = 0 (2.29) 
and, if 4 is constant, Q satisfies Laplace’s equation. 
In ferromangetic materials, equation (2.26) is not valid and it is then 


usual to treat M as a given distribution. In such materials, the funda- 
mental equations are 
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H=-—gradQ, divB=0, B= pwo(H+M). (2.30) 
Q then satisfies the Poisson equation. 
V?Q = divM. (2.31) 
The boundary conditions at a surface of discontinuity separating two 
magnetic media are (i) Q is continuous and (ii) B, is continuous. The second 
of the conditions follows from Gauss’s law by an argument analogous to 
that which led to equation (1.20). 


Problem 2.12 A long, circular, cylindrical, soft iron pipe has internal 
radius a and external radius b. Its permeability is y. It is placed with its 
axis perpendicular to a uniform field of intensity Ho. Calculate the resultant 
field. 


Solution. Taking the axis of the pipe as z-axis and x-axis in the direction 
of the applied field, we introduce cylindrical polar coordinates (p, @) as 
shown in Fig, 2.7. The potential of the uniform field is then 
—Hox = —Hg pcos $. We assume that this potential is modified by the 
induced magnetism in the pipe so that the potential in region 3 is given by 
Q = Q; = —Hopcos+Acos >/p. The added term vanishes at infinity 
and satisfies Laplace's equation as required. Inside the soft iron, we assume 
Q = 2, =(Bp+C/p)cos@, and, inside the hollow, we take 
Q = Q, = Dpcos®. Both Q, and Q, satisfy Laplace’s equation and Q, 
is finite at p = 0. Q is continuous across p = a and p = b provided 
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—Heb+4 = x Bb+E, Ba+© = Da, 


B, = Mito H, = — to OQ/dp is continuous across these two surfaces 


provided 
A 
Hota =A 5-8), (5-8) = -D. 


A, B, C, D can now be found. In particular 
jae 4utto 
GFF -G- Ie 
It follows that the field in the hollow interior of the pipe is uniform and of 
intensity 44H o/[«+ 1)? —(— 1)? a?/b?]. For large values of y, only the y? 
terms in the denominator need be retained, and the intensity approximates 
to 4Ho/[(1—a?/b)]. For large y, this will be small, showing that the 
magnetic material has a shielding effect. o 


Problem 2.13 A long circular steel cylinder of radius a is magnetized 
uniformly in a direction perpendicular to its axis, Calculate the field. 


Solution. Introduce cylindrical polar coordinates (p,#) as indicated in 


Fig. 2.8, the x-axis being parallel to M. Since M is constant, equation 
(2.31) reduces to Laplace’s equation and Q is harmonic everywhere. 


¥y 


Fig. 28 
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Inside the cylinder, we take Q = 2, = Ppcos@; outside the cylinder, 
we assume 2 = 2, = Qcos ¢/p. Q is continuous at the surface of the 
cylinder if Pa = Q/a. Inside the cylinder, B is given by equation (2.25) 
and its radial component is accordingly 
Ho(M cos @—62,/ép) = Ho(M — —P)cos ¢. Outside the cylinder, B= oH 
and the radial component is — ip 622/@p = fo Q cos $/p?. Hence, B, is 
continuous across p = a provided pio(M—P) = io Q/a*. Solving for P 
and Q, we get P = 4M and Q = 4Ma’. 

It now follows that the external field is that of a magnetic doublet of 
moment xa?M (the potential of an electric doublet of moment m is 
mcos #/2z69 p (Exercise 1, Chapter 1); for a magnetic doublet, put ¢9 = 1). 

The internal field is uniform and of intensity —$M. Note that this 
opposes the magnetization; it is termed the demagnetizing field and tends 
to destroy the magnetism. a 
2.7 Forces between Dipoles If a magnetic dipole of moment mis placed 
in a field of induction B, its potential energy is —m.B. 

Problem 2.14 Calculate the forces of interaction between two dipoles 
whose axes lie on the same plane, 

Solution. Let m, m’ be the moments of the two dipoles and let these be 
placed in relation to one another as indicated in Fig. 2.9. The forces 
acting upon m’ are equivalent to a force acting at its centre P and a couple. 
Let G be the couple and let (R, T) be the components of the force along 
the line joining the dipoles and perpendicular to this line. 
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Referring to equation (1.12), the field generated by m at m' is given by 


Ho(3m.r im 
2 B("-8) 


Thus, the potential energy of m’ in this field is W, where 


W = —m'.B= fam. mm. rm’. *) 


Homm' ro 
=F [cos(@—6’)—3 cos 0 cos 0’), 


= “oN (sin @ sin 6 ~2.c0s 0 cos 0). (2.32) 

Now suppose that m' is displaced slightly so that @ increases by d0, 0 
increases by d0’ and r increases by dr. Then P is displaced by dr in the 
direction OP and by rdé in the perpendicular direction; m’ also rotates 
through an angle d0—d6’. Thus, the work done by the forces of the system 
is Rdr+ T+d0+G(d0—d6’). This will equal the decrease in the potential 
energy —dW. But 


dW = OW 04 LL NLL 


00 00! or 
ow ow ow 
Hence Tr+G = 9" Ga R= ad 
Substitution for W from equation (2.32), now yields the results 
G= fone (sin 0 cos 0’ +2 cos @ sin 0’) 
_ 3Homm ne ay 
R= dart (sin @ sin 0’—2cos@cos @’)} (2,33) 
3utgmm' . . 
T=- ys sin(0+6'), a 


Problem 2.15 Two dipoles, of moments m and 6m, are pivoted at points 
P and Q and are free to rotate about their centres in a plane containing 
PQ. A uniform field B is parallel to this plane and perpendicular to PQ. 
Show that the position of equilibrium in which the dipoles are aligned 
with B is stable if 4nr3B > 9y9m, where r = PQ. 


Solution. Suppose the dipole m is first introduced into the field B; the 
potential energy generated is —mB cos 6@ (Fig. 2.10). The dipole 6m is then 
introduced into the combined field of m and B. Its energy in the field B is 
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Fig. 2.10 


—6mB cos ¢ and its energy in the field of m follows from equation (2.32) 
by putting m’ = 6m, 6’ = 4n—¢, and replacing 0 by 4n—0. Hence, the 
total potential energy of the system is 


3g m MA it 
w= —mB cos 0 — 6mB cos +—s- (cos 8 cos p—2sin O sin d). 


It is now easy to verify that OW/d0 = dW/dp = 0 when 0 = @ =0 
and hence that W is stationary in this configuration. Thus the dipoles are 
in equilibrium when they are Aine) with the applied field. Also, 


ew Shell eee a 

ar = mB cos 0-+—>3-(2sin 0 sin ¢ —cos 8 cos #), 
ow LA 

agF = 6mB cos p+ mn) (2 sin @ sin ¢—cos 0 cos #), 
PW — 3yym* 


2006 = Fars (Sin O sin  —2.cos 8 cos 4). 
Thus, when 0 = # = 0, 


2 
aul = mB(1—3k), oe = 3mB(2—k), 
ew ew (ew? ee 
rag a) = 3m*B*(1 + k)(2—9k), 


where k = igm/2nr°B. For stable equilibrium, W must be a minimum. 
The condition for this is that all the expressions just calculated must be 
positive, ie. —1 < k < 3. Since k is positive, the condition required is 
4nr?B > Yp19m. Oo 
Problem 2.16 The region x < 0 is occupied by material of permeability 
u. A dipole of moment m is placed at the point A(a,0,0) (a > 0) with its 
axis perpendicular to the plane face of the material. Calculate the force 
on the dipole. 
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Fig, 2.11 

Solution, The potential in the region x > 0 is assumed generated 

_ the actual dipole at A and an image dipole m’ at A’, the eptcal image oh 
in Plane face of the material (Fig. 2.11a), The potential in the region 
is assumed generated by an image dipole m" at A, regarded as 

situated in empty space (Fig, 2.11b). 

The potential is continuous across the face at a point P, provided 


m4 mar ” 
Fart ONO) + 3 080 = 75 cost —0), 

© ie. —m-+m' = —m’. (2.34) 
Dray _ The component in the direction of the normal n of the intensity at P 

_ due to m” is calculated from equation (1.12) (with & = 1) to be 

m’ 

; H, = 7-3 (8cos*0—1), 
Since B = pio H inside the material at P, 


Hom" 
B= “dans 3.00870 — 1). (235) 


Similarly, B, generated by the dipoles m, m' at P is given by 


2 B, = Ees(m-+-my3 cos?0—1). (2.36) 
fi Equations (2.35), (2.36), show that B, is continuous provided 
a m+m! = ym’. (2.37) 
Solving equations (2.34), (2.37), we find that 
{= Hole mn! = BC 
ut+l atl 
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Equations (2.33) now show that m is attracted towards m’ by a force 
3uomm! _ 3pm? w—1 
32na* —-32na* +1" oO 


EXERCISES 


1. A spherical shell having internal radius a and external radius b is 
formed from insulating material whose dielectric constant is e A point 
charge is placed at the centre. Show that the external field is the same as 
the field which would be present if the shell were removed. The charge is 
now replaced by a dipole of moment m. Show that the external field is the 
same as if the shell were absent and the dipole moment were 2m, where 


2. The planes x = 0, x = d are occupied by conductors and the region 
between them is occupied by an insulator whose dielectric constant is 
given by « = (1—x/2d)~*. The conductor x = 0 is earthed and the other 
is raised to a potential @. Show that the capacitance per unit area of the 
plate x = d is 4e)/3d. 

3. An infinite line charge q lies along the axis of an infinite circular 
cylindrical dielectric of radius a. A parallel infinite line charge —q is 
placed in the dielectric a distance f from its axis, Prove that the force per 
unit length which acts upon the first charge has magnitude 


@ [1 eis 
Qneeo| f  et+la? | 
Calculate the force per unit length on the other charge and explain why 
it is not equal and opposite to the force on the axial charge. 


4. A sphere of radius a is magnetized so that the intensity of magnetization 
at the point (x,y,z) is in the direction of the z-axis and of magnitude az, 
where « is a constant and the origin is at the centre of the sphere. Calculate 
the Poisson distribution and show that the magnetic potential outside 
the sphere is given by 2 = 2xa5P,(cos 0)/15r°. 


5. Employing spherical polar coordinates (r,0,¢), a medium has per- 
meability 41,(a/r)?, where 4, and a are constants. Show that solutions for 
the magnetic potential in the medium of the form Q = F(r)cos@ exist and 
that F(r) = Ar?+B/r. A sphere of radius a formed from this material, 
whose centre is the pole r = 0, is placed in a uniform field of intensity H. 


Show that the potential inside the sphere is given by 

‘Oa _ 3H?’ cos 

Au, + Da” 
6. Three equal magnetic dipoles are pivoted at the verti i 
e ices of an equi- 
lateral triangle, and can rotate in the plane of the triangle. Thy are 
constrained by a frictionless mechanism in such a way that their axes 
make equal angles, all in the same sense of rotation, with the inward 
bisectors of the angles at which they are situated. Show that there are four 


posto of equilibrium and that the two symmetrical positions are 


Chapter 3 
Steady Current Flows 


3.1 Kirchhoff’s Laws In this section, we consider steady electric currents 
in networks whose branches are wires meeting at points called the nodes. 
The strength of the current in a wire is measured by the rate at which 
charge flows across any cross-section of the wire. If n wires carrying 
currents i, ,i,,..-,i, meet at a node, then since, under steady conditions, 
there is no accumulation of charge at the point, 

i, +i,+...+i, = 0; (3.1) 
it is here assumed that currents directed towards the node are given a 
positive sign and currents directed away are given a negative sign. This is 
Kirchhoff’s First Law. 

Since there must be a positive component of electric field in the direction 
of the current to cause the charge to flow, the potential will decrease along 
a wire in the direction of the current. According to Ohm's Law, for any 
given wire, the potential drop along its length is proportional to the current 
flowing; thus, if i is the current and V is the potential drop, then V = Ri, 
where R is a constant for the wire called its resistance. 

Consider a closed circuit formed by branches of a network AB, BC, CD, 


...,KA. The electric potential will vary around the circuit, but must 


eventually return to its initial value upon return to the point A. It follows 
that the sum of the drops in potential along the wires of the circuit must 
be balanced by the potential rises across the terminals of the various 
batteries or other sources of current which are present in the circuit, 
It is found experimentally that the rise in potential across such a source of 
current is given by an expression E—ri, where E and r are constants for 
the source called the electromotive force (e.m.f.) and the internal resistance 
respectively, and i is the current in the sense of the potential rise. Thus, for 
any closed circuit in the network, it is necessary that 


= Ri = E(E-ri), (3.2) 

or ZDRi=ZE, (3.3) 

If it is understood that the set of resistances R includes the internal 
resistances of the sources. This it Kirchhoff’s Second Law. 

Problem 3.1 Eight equal wires, each of resistance r, form a pyramid 

ABCDE whose vertex is A and whose base is the square BCDE. Show that, 
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if current is introduced at A and removed at P, a point on BC such that 
BP = ABC, the resistance of the network between A and P is 
Tr1+A-2)/15. 


Solution. Denote the currents in AB, AC, AD, AE, DE by », w, x, y,z 
respectively (Fig. 3.1). The currents in EB, DC, BP, CP, can then be found 
(see figure) by application of the first law. Applying the second law to the 


Fig. 3.1 


circuits AEB, ADE, ADC, ABC (taking the resistances of BP and CP to 
be 4r and (1 —A)r respectively), the following equations are derived: 
2y+z—v =0; x+z—y = 0; 
2x-z-w=0; (A+ lu+dy+z—(2—Aw—(1—A)x = 0, 
Solving for v, w, x, y in terms of z, we obtain 
(1—2A)x = (2+A)z, (1 2a)y = B-Az, 
(1—22)w = (3+44)z, (1—2A)v = (7—4d)z. 
The current entering at A and leaving at P is v+w+x+ y = 1521-24), 
The potential drop along ABP is vr-+ (v+y+z)ar = 1(1+4—2)zr/(1—22). 
The result stated now follows by dividing the potential difference between 


; A and P by the total current entering the network, o 


Problem 3.2 In the ladder network shown in Fig. 3.2, the sides of every 
Square have resistance R and the ladder extends indefinitely to the right. 
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Current i enters at A and leaves at B. Find the equivalent resistance of the 
network between A and B. 


Solution. If PQRS is the rth square counting from the end AB of the 
ladder, let i, be the current in the branches PQ and RS. Then, applying the 
first law at the nodes, P, Q, we see that the currents in PS and QR are 


i 


Fig. 3.2 


i,_,—i, and i,—i,,,, respectively. Applying the second law to the circuit 
PQRSP, we obtain 


Ri, +i,-i 


red 


+i,-i,_, +i) =0, 


or ini 4, +4, = 0. 


ad | 

The characteristic equation for this linear recurrence relationship is 
2? —42+1 = 0 and the characteristic roots are accordingly 2+,/3. Its 
general solution is therefore 


i, = AQ+J/3+BQ—/3Y. 


It is clear that i, must remain bounded as r — oo and this will only be so 
if A = 0, Further, i, = i and, hence, B = i. Thus i, = i(2—,/3Y. 

Sincei, = i2-J3), it follows that the current in ABisi—i, = i(,/3—1). 
Thus, the potential drop from A to B is iR(,/3—1). This implies that the 
whole network behaves like a single resistor of resistance R(,/3—1). This 
is called the equivalent resistance of the network. o 


3.2 Heat Generation in a Network If a steady current / is flowing in a 
wire and V is the potential drop along the wire, in time t a charge it is 
transferred from the high potential end of the wire to the low potential end 
and so loses potential energy it V. This loss is compensated by a correspond- 
ing increase in the kinetic energy of the charge, which is then transferred 
by impact to the molecules of the wire and appears as heat. Thus, the rate 
of generation of heat energy is iV = i?R, if R is the wire’s resistance. 
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Problem 3.3 _A source of current has electromotive force E and internal 
resistance r. Show that the maximum rate at which heat can be generated 
in a resistor connected to this supply is E?/4r. 


Solution. If a resistance R is connected to the supply, the current i which 
flows is given by the second law E = (R+r)i. Thus, the rate of generation 
of heat in R is ?R = E?R((R+r)’. Differentiating this expression with 
respect to R and setting the result to zero, it is found that R = r makes the 
rate of heat generation a maximum. With this value of R, the rate of heat 
generation is E?/4r, as stated. o 


3.3 Steady Flows in Conducting Solids The pattern of charge flow within 
a solid conductor is specified by giving the current density vector j at all 
points. This vector is in the direction of flow and its magnitude is equal 
to di/dS, where di is the current flowing across a small surface element 
dS which is perpendicular to the flow. The rate of flow of charge across a 
surface S is |, j.m dS, where n is the unit normal to S in the sense in which 
the flow is to be calculated. If S is closed and n is outwards from the region 
T enclosed by S, this rate of flow must equal the rate at which the total 
charge inside S is decreasing (assuming there are no electrodes inside 5). 
Thus, if p is the charge density 


d 
[sma = -5 {oe 


Applying the divergence theorem to the surface integral and differentiating 
under the integral sign in the volume integral, this gives 


[ (+a i) dv = 0. 


Since I is arbitrary, this implies 
op 
a 

which is the equation of continuity. 
For a steady flow dp/dt = 0, and hence 


+div j = 0, (3.4) 


div j = 0. (3.5) 
Also, for a steady flow, Ohm’s law in the differential form 
j=ocE (3.6) 


is valid; o is called the conductivity, Finally, for a steady flow it is assumed 
that £ is derivable from a potential V by the equation 


E = —gradV. (3.7) 
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Regions from which current is supplied to a conductor, or into which 
current flows from a conductor, are called electrodes. These are usually 
assumed to have infinite conductivity and, hence, by equation (3.6), E will 
vanish within an electrode; equation (3.7) then implies that V is constant 
over an electrode. The net current supplied or received by an electrode is 
the flux of j across its surface and is called the strength of the electrode. 


Problem 3.4 A single spherical electrode of strength i is set into an 
infinite conductor of conductivity ¢. Calculate j, E and V. 


Solution. The flow will be radial everywhere and the current across a 
concentric sphere of radius r will be the flux of j across this surface, 
namely 4nr?j. Assuming steady state conditions, this must equal i, and 
hence j = i/4zr*. This leads to the equation 


t 
=r. 3.8 
j ae” G8) 
Equation (3.6) now yields 


Re eat - 39) 


and equation (3.7) then implies that 
i 


Vv =-——. (3.10) 

4xor 
Note that equations (3.8)-(3.10) have the same forms as equations (2.5)-(2.7) 
for a point charge in an infinite dielectric. Oo 


Problem3.5 A circular electrode of strength iis set in an infinite conduct- 
ing plate of conductivity ¢ and thickness t. Calculate the equations of the 
flow, 


Solution. The flow will be radially outwards from the electrode and the 
flux of j across a cylinder, coaxial with the electrode and of radius p, will 
be 2zptj. This must equal i and hence j = i/(2ntp). The equations for E 
and V are now derived as before. Thus 

i i i 
= = V=- log p, 3.11 
j 2ntp caren 2notp 4, not 8” Ge) 
where @ is the unit radial vector. o 


A similar problem is that of the infinite circular cylindrical electrode of 
strength i per unit length, set in an infinite solid conductor. The equations 
can be derived from (3.11) by setting ¢ = 1. 


5S 


Problem 3.6 Three perfectly conducting circular electrodes, each of 
small radius 4, are set into an infinite conducting plate at the vertices of 
Seah ics ts tee one Ifa current i enters at one electrode and 
currents 4i leave the other two, show that tl i resistant 
between the electrodes is vara W 
2 loe(2 
anat 8 (;). 


Solution. Since the electrodes are small, the potential at any poi 
ly point on the 
plate can be found by summation of contributions due to each (cf. point 


Fig. 3.3 


charges), Thus, if P is a point on the plate distant 
1s, Py» Pz» Ps from the three 
electrodes (Fig. 3.3), the potential at P will be given by 


i ; 
V = ~5——(log p, —}log p, —tlog p,) = aeie( 26). 
At the electrode A, p, = 5, p, = p, = a, and hence 
y= eai"*(F): 
4nat *\ 5? 
At the electrode B, p, = p, = a, Pz = 6 and, therefore, 
Vu = aeaie(). 
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3 
Thus (V,—Vpyli = aiee(3) 
This is the equivalent resistance. a 


3.4 Harmonic Functions and Images Any steady current flow is governed 
by equations (3.5}(3.7). Comparison of these equations with the earlier 
electrostatic field equations (2.2)-(2.4) reveals immediately the parallelism 
which must exist between the two theories; for, by replacing D by j, ty 
by o and setting p = 0, the electrostatic equations are transformed into 
those for a current flow. However, the boundary conditions for flow 
problems take a distinctive form and these will now be stated: (i) over the 
surface of an electrode, V is constant and, either V takes a given value 
(electrode’s potential known), or the flux of j out of the surface takes a 
given value (electrode’s strength known); (ii) over an insulated boundary, 
Jj, = 0, since there is no flow from the conductor into the insulator; also, 
at a surface separating two conductors having different conductivities 
@, and ,, V and j, are usually assumed to be continuous; j, must be 
continuous if there is to be no accumulation of charge at the interface, 
but the continuity of V depends upon there being no chemical interaction 
between the two conductors, which might lead to a contact potential 
difference being generated. 

Eliminating j and E between equations (3.5)-(3.7), it is found that V 
satisfies the equation 


div(o grad V) = 0. (3.12) 


Thus, over any region for which g is constant, V satisfies Laplace’s equation 
and the methods of harmonic functions and images previously described 
become applicable. 


Problem 3.7 A conducting sphere of radius a and conductivity o, is 
embedded in an infinite conducting medium of conductivity o,, and 
modifies a steady current flow which would otherwise have a uniform 
current density j,. Show that the current density in the sphere is 
30, jo, +2c,) and calculate the density of charge on the spherical 
interface. 

Solution. This problem is mathematically similar to the electrostatic 
problem of a dielectric sphere placed in a uniform field (Problem 2.6), 
Taking the pole of spherical polar coordinates at the centre of the sphere 
and the z-axis in the direction of j,, the potential associated with the 
uniform flow will be V = —j,rcos 0/c, (for this corresponds to a uniform 
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field of intensity E = j,/o, and, hence, j = o, E = j,). We now assume 
that the potentials inside and outside the sphere are given by 
= = es a 
V, = Broos@, =v, = Pa 
respectively. The constants A, B are now chosen to ensure that 
J, = —(@V/dr) and V are continuous across the surface r = a; as in 
Problem 2.6, we find that 


A= 212% Tio Bie 
o,+2¢,° 0," o, +20," 


Thus, we have shown that the field within the sphere is uniform and of 
intensity E = 3j,/(c, +2c,); the associated flow is accordingly also 
uniform and of density j = o, E = 36; igo, +20,). 

Calculating the radial component of E just inside and just outside the 
sphere, we find that 

E,, = —0V,/ér = 3j, cos 8a, +20,), 
E,, = —@0V,/dr = 30, jg cos 0/o,(a, +20). 
Assuming that the dielectric constants of the conductors are unity, equa- 


tion (1.20) shows that the surface density of the charge accumulating on 
the interface is 


3(¢,-93) 
8(E,,—E,,) = ale, $2a,) e008 8. 
(We are assuming the equation div D = p remains valid in these circum- 
stances; see § 5.6.) Oo 


All the electrostatic problems solved in Chapters 1 and 2 can be restated 
as current flow problems and their solutions remain unaltered in form. 
The only essentially new feature by which current flow problems differ 
from those relating to electrostatic fields in dielectrics, is the boundary 
condition j, = 0 over the insulated surface of a conductor. We shall 
accordingly close this section by giving two examples of problems in 
which this condition arises and which are solved by introducing image 
electrodes. 


Problem 3.8 A semi-infinite solid of conductivity ¢ is bounded by a 
plane face. A pair of small spherical electrodes of radius 6 are set into the 
solid a distance 2a apart and at the same distance a from the plane face. 
Calculate the equivalent resistance between the electrodes. 

Solution. Suppose that a current i enters the conductor via electrode A 
and leaves via electrode B (Fig. 3.4). Imagining the whole of space to be 
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ent 


filled with the conducting medium, we introduce image electrodes of 
strengths +iand —iat A’ and B' respectively, A’ being the optical image of 
A in the boundary plane and B’ being the optical image of B. At any point 
P on the boundary plane it follows from the symmetry of the electrode 
system with respect to this plane that the resultant current density vector 
of the flow generated by the four electrodes is parallel to the plane. Hence, 
if m is the unit normal to the plane, j.m = 0 and the boundary condition 
is satisfied. We conclude that the flow generated by the four electrodes is 
identical with the actual flow in the semi-infinite region occupied by the 
conductor. 


> 


y 


4, 


Fig. 3.4 


Employing equation (3.10) the potentials at the surfaces of the electrodes 
A and B, can now be found; they are given by 


sofa odd [5-7 
"= gel jes —ayae | - ae 3 ZJua)’ 
i be deteihve dd ee | 
= ge Fatty il -5tam [ 
The equivalent resistance between the electrodes is therefore 
V,-V, 1] 1 
aa 1 aaa Q 
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Problem 3.9 Electrodes of small radius 5 are set into a conducting 
circular disc of radius a, thickness t and conductivity ¢. The electrodes 
lie on the same diameter at equal distances f from the centre. Find the 
equivalent resistance between them. 


Solution. Suppose a current i enters the disc at A and leaves at B (Fig. 
3.5). Imagine the disc material is extended to infinity and that additional 
electrodes of strengths i and —i are introduced at the points A’, B’, which 


Fig. 3.5 


are the inverses of A, B in the circle, The resultant current density at a point 
P on the edge of the disc due to the two electrodes at A and A’ is found by 
the use of equation (3.11) to have a component in the direction of the 
radius OP equal to 
i i 
Intro Fey 
where r= AP, r’ = A‘P. Putting r’ = ar/f, cos = (a? +r? —f?)/2ar, 
cos @ = (a? —r? — f*)/2rf, in this expression, it will be found to reduce to 
i/2nta. An exactly similar argument shows that the same component of 
the current density generated by the electrodes at B, B’ has value —i/2ta. 
Thus, the flow generated by all four electrodes satisfies the boundary con- 
dition j, = 0 at the edge of the disc and accordingly represents the actual 
flow in the disc. 
The potential of the electrode A is now calculated (equation (3.11)) to be 


i= Zeal —Iee+1082/-W0e(%—1) +100(*+1)] 
< imaloe| si 


~ 2xot ia? —f?) |" 


cos 0, 


The potential V, of the electrode B is the negative of this. The equivalent 
resistance between the electrodes is accordingly 

as! 2f(a* +f?) 
7 = Sloe AS : o 


EXERCISES 


1. Conducting wires, all of equal resistance r, form the edges of a cube. 
Prove that the equivalent resistance between two adjacent corners is 
Tr/12. 


2. Five small circular electrodes, each of radius 6, are set into an infinite 
plane sheet of conductivity a and thickness ¢, at the vertices and centre ofa 
square of side a. A steady current 4i enters the sheet at each of the vertices 
and a current i leaves at the centre. Show that the equivalent resistance 
between the four vertices and the centre is 


ae ne 

Brot 8 \ 205)" 
3. A perfectly conducting sphere of radius a is surrounded by a concen- 
tric spherical shell of inner radius a and thickness Aa, and of conductivity ¢; 
the rest of space is filled by a medium of conductivity a’. The conducting 


sphere is maintained at a potential ¢ above that of points at an infinite 
distance. Show that the current flowing from the sphere is 


4naaa'(A+1)o/(a +o’). 


4. A uniform conducting medium extends to infinity in all directions 
and has conductivity ¢,. It supports a uniform flow whose density is j,. 
A spherical cavity of radius a is cut in the medium and is filled with 
material whose conductivity at distance r from the centre is ¢, a/r. Assum- 
ing that the potential within the sphere is of the form V = y(r)cos0, 
where r, @ are spherical polar coordinates, show that the ratio of the total 
current through the sphere to the current which flowed through the same 
region before the cavity was cut, is 3/(\/2+ 1). 


5. Two small circular electrodes, each of radius 6, are a distance 2f 
apart in an infinite plane sheet of thickness t and conductivity o. A circular 
hole of radius a(< f), whose centre is midway between the electrodes, is 
punched in the sheet. A current i enters at one electrode A and leaves at 
the other B. By placing image electrodes of strengths +i and —i at the 
points inverse to A and B in the circle, show that the effective resistance 
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between the electrodes is 
1 2? +a?) 
zat Oa ¥ 


6. Current flows between two perfectly conducting circular electrodes of 
radius a in an infinite plane sheet of conductivity ¢ and thickness 1. 
Iftheir centresarea distance b(> 2a) apart, show that the resistance between 


them is 
not °°| 2a \4a?— bi 


(Hint: Refer to Problem 1.22.) 
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Chapter 4 
Steady Electromagnetic Fields 


4.1 The Biot-Savart Law In this chapter, we shall be concerned with 
the magnetic field associated with steady currents of the types studied in 
Chapter 3. The electric field responsible for the flow of charge together 
with the magnetic field are regarded as two aspects of a single unified field 
called an electromagnetic field. However, in the steady conditions assumed 
in this chapter, no interaction between the two component fields takes place 
and they may, accordingly, be studied separately. 

The theory will be based upon the Biot-Savart Law: Let ds be any 
vector element of a wire carrying a current i (ds is taken in the same sense 
as i) and let r be the position vector, relative to the element as origin, of 
any point P (Fig. 4.1). Then the contribution of the charge flowing in the 
element to the magnetic intensity at P is given by 

idsxr 
dH = (4.1) 


This means that dH has magnitude ids sin 0/4nr?, where 0 is the angle 
between ds and r, and its direction is that in which P would move if it 
were rotated in the positive sense (right-hand screw) about an axis along 
ds. 


Problem 4.1 Calculate the magnetic intensity at any point on the axis 
of a circular wire of radius a carrying a current i. 
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Solution. P is a point on the axis of the wire distant x from its centre O 
(Fig. 4.2). The contribution to the magnetic intensity at P of the element 
QQ’ of the wire has magnitude i ds/4nr?, where ds = QQ and QP =r: 
the direction of this contribution is perpendicular to both QQ’ and QP, 
ie. at an angle 3x ~« to OP. The axial symmetry indicates that the resultant 
field at P will be directed along OP and it follows that only the component 
idssin «/4zr? of the contribution of QQ’ need be retained. Integrating 


Fig. 4.2 


these contributions around the circle, we obtain the net intensity at 
Pas 


. isin a iasina ia? 
A= ag |= F- =3aae 02 
ff In particular, at the centre O the field has intensity i/2a. o 


_ Problem4.2 A steady current i is flowing around a square circuit of 
side 2a. Prove that the magnetic intensity at the centre of the square has 
_ magnitude ,/2i/na. 


Solution. Let QQ’ be an element dx of the side AB of the square distant x 
from its midpoint (Fig. 4.3). The current i flowing in this element contri- 
butes to the field at the centre O an intensity idx cos 0/(4nr?) in a direction 
perpendicular to the plane of the circuit. Putting x = atan0, dx = 

asec? @d0, r = asec, this contribution takes the form icos 0d0/(4za). 


By integration over the range —4n < 0 < 4n the net contribution of the 
side AB is found to be 


i 


Fig. 4.3 


The remaining three sides each make the same contribution, so that the 
result follows. ‘ ; o 

By extending the range of the above integration to ud <0<h, 
the intensity at a distance a from a very long straight wire carrying a 
current i is found to be given by 


ha, (4.3) 


Problem 4,3 Wire is wound closely in a spiral fashion to form a circular 
disc. Current is sent through the wire so that i is the current across unit 
length of any radius. Show that the field intensity at any point on the axis 
of the disc is 4i[cosh™ ‘(sec #)—sin«], where « is the angle subtended at 
the point by any radius of the disc. 

Solution. The current flowing around a ring-shaped element of the disc 
(Fig. 4.4) of radius r and width dr is i dr, It follows Roc easter ae 
the intensity contributed by this element at P on the axis is tir’ arr +x773, 
Putting r = xtan@, dr = xsec?0d0, this reduces to titan?0cos0d0. 
Hence, the net intensity at P due to the whole disc is 


4i [ tan? cos @.d0 = i , (sec @—cos 0) d0 
= tiflog(sec a + tan )—sin a]. 
This is equivalent to the result stated. 
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It can be proved that the law (4.1) is equivalent to the statement that 
the field generated by a circuit C carrying a current i can be derived from 
magnetic potential 2 = im/4n, where « is the solid angle subtended at a 
point of the field by C. Thus 


H = —grad (2). (4.4) 


It should be noted that this potential is not single-valued, since the 
solid angle is arbitrary to the extent of an added multiple of 4x, Thus, 


P 


A 


| 


= 


Fig. 4.4 


if P is taken around a closed path linking C, the continuous variation of w 
will cause its value to change by 4z. The implication is that the field is 
not conservative; the fact is that such fields can be made to perform work, 
the energy necessary being supplied by the current source. The operation 
of the electric motor is based upon this principle. 


Problem 44 Calculate the magnetic potential due to a long straight 
wire carrying a current i, 


Solution. The potential Q can be calculated by use of the result (4.4), 
but it is more easily derived from equation (4.3). 
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Since the field is two-dimensional, @ will depend on the cylindrical 
polar coordinates (p, #) alone (the wire is being taken as z-axis). The polar 
components of the field intensity are given by 


4 1 
ODY ee pop 2np° 
It now follows that Q is independent of p and that 
ip 
=> (4.5) 


Since ¢ is multi-valued, we again note that @ is multi-valued also. Also, 
since the reference line @ = 0 can be chosen arbitrarily, this potential is 
indeterminate in respect of an added constant. Qa 


4.2 Force on a Current If a wire carrying a current i is placed in a 
magnetic field of induction B, each element ds of the wire experiences 
a force dF given by 

dF = idsxB. (4.6) 
Note that the force is perpendicular to both the current and the field. This 
formula remains valid if the current is embedded in magnetic material. 

A particular case is the force between a pair of long straight parallel 
wires carrying currents i and i’ and distance a apart. Using equation (4.3), 
the force is found to be y4, ii’/2na per unit length and is one of attraction 
if the currents are in the same sense and of repulsion otherwise. 


Problem 4.5 A current i is flowing around a circle of radius a. A dipole 
of moment m is placed on and along the axis of the circle at a point 
distant c from the centre. Show that the force acting on the circuit is 
34, ima? cla? +c7)~*, 


Solution. Let ds be an element of the circuit at a point P (Fig. 4.5) and 
let r be the position vector of P with respect to the dipole at O. The field 
at P due to the dipole follows from equation (1.12); substituting in equation 
(4.6), it is found that 


i[ 3m.r dsxm 
dF = fall 2 asx 5. 
ds x m is directed along the radius PC and the corresponding component 
of dF is cancelled by a similar component acting on the element of the 
circuit diametrically opposite to ds. ds x r is directed along the the perpen- 
dicular to OP lying in the plane OPC and its contribution to dF can be 
split into two components (i) rds cos @ along PC and (ii) rds sin @ parallel 
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to OC. The first of these components is again cancelled by a similar 
component arising from the element opposite P. The second component 
can be integrated over all elements of the circuit to give a resultant force 
— | Hol 3mcosé, 
F= ae ae dssin0 
along OC. Since r, 0 do not vary as P describes the circle, this gives 
3g ima sin 8 cos 8/2r°, which is equivalent to the result stated. Oo 


Fig. 4.5 


Problem 4.6 _Find the force of interaction between a current i flowing 
ina circular wire of radius a and a current i flowing in a long straight wire 
which lies in the plane of the other wire and is distant (> a) from its 
centre. 


Solution. Consider the element PP’ of the circular wire subtending an 
angle dO at its centre O (Fig. 4.6). The magnetic induction at P due to the 
current i’ has magnitude yi, i'/2nPN = Hy t'/2n(c—acos 8) (equation (4.3)) 
and its direction is perpendicular to the plane of the two wires. The force 
exerted upon PP’ is therefore Hg ii'a d0/2x(c —acos 6) in the direction OP. 
Because of the symmetry about OA, only the component of this force 
along OA contributes to the resultant force acting upon the circle and this 
has magnitude , ii’a cos 0 d0/2x(c —acos 6). Integrating this over the 
range 0 < 0 < 2r, the resultant force of attraction of the circle towards 
the straight wire is found to be 

m ce 
Hy it [zea] 
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(Note: The integration can be performed by expressing the integral as a 
contour integral around the circle |z| = 1 in the complex plane, or by 
changing the variable of integration to r = tan 40.) oO 


Fig. 4.6 


4.3 Potential Energy of a Current Loop If a steady current i is caused to 
flow in a loop of wire and this is placed in a steady magnetic field whose 
induction (in the absence of the loop) is B, then the potential energy of the 
system is —iN, where N is the flux of B through the loop (more precisely, 
the flux of B through any surface S bounded by the loop). This statement 
remains valid if the current is embedded in magnetic material. 

If the wire is displaced and the current i is maintained constant, the work 

done by the forces acting upon the wire will equal the decrease in the 
potential energy. 
Problem 4.7 If the current i lies in a plane and encloses an area A and if 
the field B is uniform, show that the system of forces acting on the loop is 
equivalent to a couple of moment iAa x B, where a is a unit normal to the 
plane of the loop. 


Solution. Suppose the system of forces acting on the loop is equivalent 
to a force F acting at some point O in its plane and a couple G. If the loop 
is translated (without rotation) through a small displacement da, the 
work done by the force system is F.da. But the flux N of B through the 
loop remains unchanged and, hence, F.da = 0. Since this equation is 
valid for arbitrary da, F vanishes and the force system is equivalent to a 
couple. 
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~ Next suppose the loop is given a small rotation d@ about O. Then the 
normal will undergo this rotation and will change to m +d@ xa. Thus, 
the flux of B through the loop will change from AB.n to AB.(n+d0 xn), 
ie. will increase by AB. d@ xn. This means that the potential energy will 
decrease by iAB.d0x n= iAnxB.d0. But the work done by G is 
G.d0 and, hence, G.d@ = iAn x B.d0. This being true for arbitrary d0, 
we conclude that G = iAnx B. 

It should be noted that this is exactly the couple which the field B would 
exert upon a magnetic dipole of moment iAa. This is an important result; 
namely, any small current loop i enclosing a plane area dS whose unit 
normal is # behaves like a magnetic dipole of moment idSn. Qa 


4.4 Method of Images If a steady electromagnetic field is generated by 
electric currents, then over any region devoid of currents, the equations 
(2.28) or (2.30) which governa field due to permanent magnets are applicable. 
Unlike the magnetostatic case, however, the potential Q will be multi- 
valued. In particular, within a material of constant permeability, Q will 
satisfy Laplace's equation and the method of harmonic functions and 
images are applicable. 

Problem 4.8 Soft iron of permeability 1 occupies the whole of space to 


one side of a plane. A current i flows in a long straight wire which is outside 
the iron and parallel to its plane face. Show that the wire is attracted 


towards the iron by a force j1, i?(u—1)/4za(u+ 1) per unit length, a being 
the distance of the wire from the iron. 


Solution. The method is similar to that employed in Problems 2.8 and 
2.16. Outside the iron, the field intensity is assumed to be generated by the 
actual current i through A and a parallel current i through A’, the optical 
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image of A in the plane face (Fig. 4.7a). Inside the iron, the field intensity is 
supposed generated by a current i” in the actual wire through A (Fig. 4.7b). 
Then, the potential of the field outside the iron at P on the plane face is, 
by equation (4.5), [—i(x—6)—i'0]/2n. The potential at P for the field 
inside the iron is —i”(x—6)/2n. Since an arbitrary constant can be added 
to these functions, the condition that the potential should be continuous 
across the plane face is 
i-i =i". (4.7) 
Using equations (4.3) and B = y1, H, it is found that the component of 
B at P normal to the plane face generated by i and 7’ is given by 
B, = —Ho(i+i')sin 0/2xp, Inside the iron at P, B, = —HHy i” sin 0/2xp. 
Hence, B, is continuous across the face provided 5 
iti = pi”. (4.8) 
Solving equations (4.7)and (4.8) we get i’ = (u—1)i/(u+1),i" = 2i/(u+ 1). 
The force acting upon i is that due to the image current i’ and, using the 
result derived on page 67, this is now found to be the result stated. 


Problem 4.9 A long circular cylindrical cavity of radius a is drilled — + / 


through a large block of iron of permeability 4. A current i flows ina 
straight wire parallel to the axis of the cavity and returns along a second 
similar wire. The wires are on opposite sides of, and at the same distance f 
from, the cavity axis. Calculate the force between them. 


Solution. A right-section of the cavity is shown in Fig. 4.8, A and B 
being the points in which the plane of section is intersected by the two 
wires. We introduce parallel image currents +/’ to meet the plane in the 
points A’, B’ inverse to A, B in the circular boundary of the cavity, and 
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assume that the magnetic intensity in the cavity can be generated by the 
system of four currents. The intensity in the iron we shall assume can be 
generated by currents +i” at A, B. 

Then, the normal component H, along PO of intensity at P in the cavity 
generated by the currents at A, A’ is (isin ¢/2mr)+(i’sin 0/2zr’). Since, by 
the sine rule, sin ¢/r = sin /r’, this component reduces to (i+i’)sin0/2nr’. 
Thesame component of Hin the iron at P generated by i” at A isi” sin 0/271’. 
The corresponding contributions of the currents at B, B’ will be 
—(i+/)sin 0,/2nr, and —i’sin@,/2nr,, where 0, = angle PBO and 
1, = BYP. It now follows that B, is continuous at P provided 

i+? = pi", (4.9) 

The magnetic potential at P in the cavity due to i, i at A, A’ is 
—[i0+i(—9)]/2x = —[i0+i(n-0+Y]/2n. Similarly, the potential 
contributed by —i, —/ at B, B’is (10, +i',)/2n = [i0, +? —0,)|/2n.The 
net potential due to the four currents is accordingly (i— i')(0, — 0)/2n (neg- 
lecting an added constant). Similarly, the potential at P in the iron due to 
the currents +i” at A, B is calculated to be i"(0, —0)/2x. The condition for 
the potential to be continuous at P is therefore 


i-f =i". (4.10) 
Solving equations (4.9), (4.10), we find that i! = (u—1)i(u+1), 
i” = 2i/(u+1), It now follows that the resultant force acting on unit length 
of the current iat A due to the currents at A’, —iat Band —i’ at B’, which 
is directed along OA, is (see page 67) 
mela Ph ghoh fo) iE iL ak 2a 
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4.5 Ampére’s Law If the steady flow of charge is not confined to wires 
but takes place in a three-dimensional conducting medium and is specified 
by a current density j, the intensity H of the magnetic field generated is 
given by the equation 

curl H = j. (4.11) 
This equation follows from equation (4.4) thus: if C is a closed path 
embracing the current i and it is described in a positive sense with respect 
to i, equation (4.4) implies that 


i 
f aa = -7z,(¢], 
c 
where [] is the increment in @ for one circuit of C. But [w] = —4z and 
hence 
n2 


ped =i. 
z 


This is Ampére’s Law. Assuming this law to be valid for the steady flow j, 
since the current through C will be given by the flux of j across a surface $ 
bounded by C, we have 
$H.dr = | j.nds. (4.13) 
ce s 
Application of Stokes’ Theorem (L. Marder, Vector Fields, p. 50) now 
leads to equation (4.11). 
To summarize, the equations of the magnetic field in a medium of 
permeability due to a given current distribution j are 


curl H = j, div B = 0, B= py. (4.14) 
In regions where j is zero, a potential Q exists and 
H = —gradQ. (4.15) 


Problem 4.10 An infinitely long straight rod has a circular cross-section 
of radius a. A current i flows in the rod, its distribution over any cross- 
section being uniform. Calculate the magnetic field generated. 


Solution. Taking axes Oxyz, with Oz along the axis of the rod, and intro- 
ducing cylindrical polar coordinates (p, ¢, 2) in the usual way, let(B,, B,, B.) 
be the corresponding radial, transverse and z-components of B. The axial 
symmetry requires that all these components shall be functions of p alone, 
Hence, constructing a right circular cylinder, coaxial with the rod, of 
radius p and unit length, thé fluxes of B out of the two circular ends will 
cancel; also, the flux of B out of the curved surface will be 2xpB,, But the 
flux of B out of any closed surface always vanishes and, hence, B, = 0. 
Thus H, = 0. 

The current through one of the circular ends of the cylinder is i if 
p > aand is ip?/a’ if p < a. Thus, if H, is the component of H tangential 
to the circumference of this end, then H, is constant and Ampére’s law 
applied to this circle gives 

2npH, = ip*/a, p <a, 
=i, p>a 
Hence, inside the rod H, = ip/2na?, and outside, H, p= W/2np. 

It follows from the previous argument that H,, H,, H, are functions of 
x, y alone. Hence, the x- and y-components of curl H are dH./éy and 
—6H,/Ox respectively, But, equation (4.11) requires that these components 
should vanish and, therefore, H_ is constant, Since H must be zero at an 
infinite distance, this means that H, = 0. 
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We have proved, therefore, that the field outside the rod is identical with 
that due to a current i flowing along its axis. Oo 


Problem 4.11 A current of constant density flows in an infinite, plane, 
thin conducting sheet, the lines of flow being parallel. Calculate the field 
outside the sheet. 


Solution. Take axes Oxyz so that the x-axis lies along a line of flow and 
__ the z-axis is normal to the sheet. Then H can only depend upon z and the 
components of curl H are therefore (—0H,/0z,0H,/0z,0). Since curl H 
vanishes outside the sheet, H, and H, must both be constant. Hence, 
div B = 0B./dz = 0 and, thus, B. is also constant. In the region z < 0, 
suppose B. = B,; then, in the region z > 0, B, = —B,. But the normal 
component B. is continuous across the sheet and it follows that B, = 0. 
Thus H, = 0. 

Now construct a rectangular circuit in the yz-plane whose sides are the 
lines y = +y,z = +2, and apply Ampére’s law. If is the current density 
in the sheet and t is its thickness, the current flowing through the circuit 
is 2jtyy, The tangential line integral of H around the circuit is 2y,(H, yay) 
where H, = H,, for z <Q and H, = H,, for z > 0. Ampére’s law shows 
that H,,—H,, = jt, But H, must have the same magnitude in the two 
regions and, hence, H,, = 3jt and H,, = —3jt. 

Byasimilar application of the law to a circuit in the xz-plane whose sides 
are the lines x = +x,, z = +2, we prove that H,, = H,, = 0 (there 
is no current flow through this rectangle), 

We have shown, therefore, that the magnetic field in the region z < 0 
is uniform, of magnitude 4j¢ and parallel to Oy. The field in the region 
z > Ois the same, except that it is in the reverse sense. o 


EXERCISES 


1, A current circulates around a wire bent into the shape of a semi- 
circle with its diameter. Show that at a point on a line through the centre 
perpendicular to the plane of the semi-circle, the intensity makes an angle 
@ with the line given by tang = (2tana)/n, where 2x is the angle sub- 
tended by the diameter at the point. 


2. A wire carrying a current i is in the shape of the plane curve having 
polar equation r = f(0). Show that the magnetic intensity at the pole 
has magnitude 

i fa 

4njor° 
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4 . v i, ~ ‘ . 
Deduce that the field at the focus of an elliptical current has intensity 
i/l. where | is the latus rectum. 
3. Accurrent is flowing in a square circuit having side of length a. A side 
of the square subtends an angle 2x at a point P on the perpendicular to the 
plane of the square through its centre. Show that the magnetic intensity 
at P has magnitude (4isin « tan?a)/na. 


4. A current i flows around a circle of radius a and another current i! 
flows in a long straight wire which is perpendicular to the plane of the 
circle and a perpendicular distance c from its centre. Show that there is a 
couple tending to bring the two wires into the same plane, of moment 
Mp i'a?/2c oF Hy ii’c/2, according as ¢ > a or ¢ < a. 
5. A current i flowing in an infinite, straight wire is at a distance S from 
and parallel to the axis of an infinite circular cylinder of iron of permea- 
bility . The radius of the cylinder is a(</f). Show that the presence of the 
iron reduces the magnetic intensity over the region it occupies by a factor 
2/(u-+ 1) and that the wire is attracted to the iron with a force 
Hoa? ?(u~1) 

2nf (f? —a*\u+1) 

per unit length. 
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Chapter 5 
General Electromagnetic Field 


5.1 Electromagnetic Induction If a loop of wire C is placed in a varying 
magnetic field, an e.m.f. is induced in the wire and causes a current to flow 
around it. The e.m-f. V is proportional to the rate of change of the flux of 
magnetic induction N through the loop; if SI units are employed, we can 
write 
dN 

V= a (5.1) 
aresult which is called Faraday’s Law. The negative sign is correct provided 
the positive senses of the e.m.f. around the circuit and of the flux through 
the circuit are related by a right-hand screw rule; the existence of this sign 
is frequently called to attention by reference to a back emf. 

If S is any surface bounded by C, and nis its unit normal, then 


Ne J B.nds, (5.2) 


Since div B = 0, N is independent of the form of S, 

Thus, ifa plane loop of wire enclosing an area A is rotated with angular 
velocity win a uniform magnetic field B perpendicular to the axis ofrotation, 
the flux of induction through the loop when the normal to the plane of the 
loop makes an angle wt with Bis N = ABcos wt, Substituting in equation 
(5.1), the e.m.f. induced is calculated to be V = ABwsinwt. An emf. 
varying with ¢ in this manner is called an alternating emf. 


Problem 5.1 A magnetic dipole of moment m is placed with its axis 
along the axis of a circular wire of radius a at a distance c from its centre. 
Ttis then moved along the axis to this centre. If R is the resistance of the wire 
and its inductance is negligible, show that a charge 


wae fee Vs 
2Ra| \a? +c? 


Solution. At time t during the motion of the dipole, let N be the flux 
of B due to the magnet through the wire. Then dN/dt is the e.m.f. induced 
in the wire and, by Ohm’s law, the current which flows is i = daN/Radt. 
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is caused to circulate. 


Thus, the charge which circulates during a time interval (t,,t,) is equal to 
es T (tan 1 N.—N 
q= [ idt= al a = R '2—N;). 


Let P be a point in the plane of the wire having spherical polar coordin- 
ates (r, 0) relative to the initial position of the magnet as pole (Fig, 5.1). 


Fig. 5.1 


It follows from equations (1.11) (replacing 1/e, by 1.) that the component 
of Bat P normal to the plane of the wire is 1, m(3 cos? @— 1)/(4nr*). Thus, 
if CP = y, the initial flux of B through the wire is 


fe ta bi 29 , 
ae ep eee 1), 2zy dy. 


Putting y = ctan0, dy = csec?0d0, r = csec®, this flux is found to be 
N, = 4,9 a?m(a? +c?)~?. The final flux is found by putting c = 0 to be 
N, = Hy m/2a. The value of q now follows. Qo 


5.2 Inductance Between Circuits If the loop C of the last section is 
situated in a magnetic field generated by a current i’ in another loop C’, 
and we assume that the magnetization induced in any materials present 
is everywhere proportional to H (equation (2.26)), then the induction B 
at a given point due to the current i’ will be proportional to i’. Thus we can 
write N = Mi’. M is called the mutual inductance of C and C’. It may be 
proved that this relationship between C and C’ is a reciprocal one and 
that, if both circuits are immersed in an infinite medium of permeability 1, 


hee a (53) 
4u r 
12 


where ds, , ds, are vector elements of the circuits and r is the distance 
between them. 
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Problem 5.2 Calculate the mutual inductance of the two circuits shown 
in Fig. 4.6. 


Solution. Taking OA as x-axis and the perpendicular to this line through 

CO in the direction of the current i’ as y-axis, the induction due to this 

current at a point (x, y) is of magnitude B = yi'/2n(c—x). Bis everywhere 

perpendicular to the plane of the circular circuit. Integrating B over the 

area enclosed by the circle, we obtain for the flux of induction through 
this circuit, 

1 fa ha? = x2) 
N= tet , a | aa = 


Is =x 


It can be verified by differentiation that 
ox @ 


a(c—x) 


Mea 2 ae = esin-!=— J/(c?—a?)sin —J(a@—-x?). 
Thus, N= Mg i[e—J(c?—a*)] and the mutual inductance is 
Hole—V(c? —a*)}. o 


Problem 5.3 Two equal, circular loops of radius a, lie opposite each 
other, a distance c apart. Show that their mutual inductance is 


an cos 0 dO 


M=4a,) =>. 
27 Ho lo (c? +2a* —2a*cos 0)* 


G 
Ifc is large, deduce that M = 1 na*/2c’. 


Solution. Elements ds, , ds, of the two circuits are shown in Fig. 5.2 as 
PP’, QQ’ respectively, Taking their scalar product, we find 


wc) 


ds, .ds, = a’cos(0—) d0d@. Relative to the axes Oxyz, the coordinates 
of P are (acos 0, asin 6,0) and of Q are (acos #, asin ¢,c). Hence 


PQ =r = [c?+2a?—2a’cos(—¢)]*. 
Equation (5.3) now leads to the result 


eds 2m (2x a?cos(0—4) 
~ 4 jo Jo [c? +2a7—2a?cos(0— A)? 


0 do. 


But 


a cos(@—) dO acy cos 0 dO 
jo [e?-+2a*—2aPcos(— A]? J, (C7 ++ 2a? = 2a cos 6) 


This integral is independent of #; the integration with respect to @ accord- 
ingly simply introduces a factor 2x and the result stated then follows, 

If a/e is small, (c? +2a?—2a?cos 0)~* = c~'[1—a*(1—cos)/c?] to 
order (a/c)?. The integral upon which M depends has therefore the 
approximate value za*/c?. Hence, M = pt na*/2c?. Oo 
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5.3 Circuits with Discrete Elements In this section, we shall study electric 
circuits in which elements called resistors, capacitors and inductors are 
joined together by wires having negligible resistance. 

Resistors and capacitors have been sufficiently described in earlier 
chapters. An inductor is usually a coil of wire wound on a cylindrical 
former of magnetic material. Suppose a current i is flowing through the 
coil; assuming that the magnetic induction generated is everywhere 
proportional to i, the flux of B through each turn of the coil will also be 
proportional to i. Thus, the total flux of induction through the coil can be 
taken to be Li, where L is a constant for the element called its self-induct- 
ance. By Faraday’s law, if i varies an e.m.f. V will be induced in the coil 
given by the equation. 


Vile (5.4) 


The negative sign indicates that this e.m.f. acts so as to oppose the change 
in the current. 

If the resistance of an inductor cannot be neglected, this will be repre- 
sented in a circuit diagram by a separate resistor in series with the coil. 


Problem 5.4 An inductor having self-inductance L and resistance R is 
connected in parallel with a capacitor C in series with a resistor R’. If the 
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combination is connected to a source of constant e.mf. E, find the condi- 
tions satisfied by L, C, R, R’ if the current taken is constant. 


Solution. The circuit is shown in Fig. 5.3. Let i, , i, be the currents flowing 
in the two branches at time ¢ after connecting the supply and let g, —q be 
the charges on the plates of the capacitor at this instant. Then, the drops 


i 


Fig. 5.3 


in potential between the ends of the branches are Ri, +Ldi,/dt and 
Ri, +q/C. These must both equal the e.m-f. applied by the supply; hence, 
since dq/dt = i, , the following equations are valid: 
ey ee See ae 
E=Ri,+lL = RT 4+. (5.5) 
Solving for i, and q under the initial conditions i, =0,q=Oatt=0, 
we find that 


sala) enfin) 0 


Hence i, = dq/dt = (E/R')exp(—t/CR’). 

If i = i, +i, is to be constant, the two exponential terms must cancel. 
The conditions for this are R = R’ and R/L = 1/CR’, ie. L = CR?. The 
constant current flowing is then E/R. o 


Problem 5.5 A pair of identical circuits, each containing a capacitor C 
and a coil of inductance L (negligible resistance) in series, are so arranged 
that the coils are magnetically coupled and have mutual inductance M. 
At t = 0, the currents in both circuits are zero and the capacitors have 
charges Q and 0. Determine the charges at time t > 0. (Assume M < L.) 
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Solution. Employing the notation of Fig. 5.4, and equating the potential 
drops around the circuits in the senses of the currents to the e.m.f.s in th 
circuits, we get the equations 

4 = -LDi,-MDi,, 4 = -LDi,—MDi,. 

Since i, = Dg, ,i, = Dq,, these equations are equivalent to 
(LD? +1/C)q, + MD"q, = 0, 
MD*q,+(LD? +1/C)q, = 0. 

Eliminating g,, we arrive at the following equation for qj: 
[((L—M)D? + Vai [(L+M)D?+1/C]q, = 0. This equation has character- — 
istic roots +jA, +ju, where 2? = 1/C(L—M), »? = 1/C(L+M), and the s 
general solution for q, can accordingly be written, 

q, = PcosAt+P’sin At+ R cos yt + R’sin pt. 
Substituting this in the first of equations (5.7) and integrating twice, we 
find that g, = —P cos At—P’sin At + R cos pt + R’sin pt. 


(5.7) 


Fig. 5.4 


The initial conditions are t = 0, g, = Q, gq, = 0, Da, = Dg, =0. © 
To satisfy these, it is necessary that P = R = 40, P’ = R’ = 0; thus, 
4, = 4Q(cos pt-+cos 21), g, = 4O(cos ut—cos At). This is the solution 4 
required. It can also be written in the form . 

4, = Qcos(A—pltcosh(A+p)t, q, = Qsin}A—p)rsin}(A+ we. 
Then, if M is small by comparison with L, A and y differ by a small quantity 


only and both q,, 4, execute oscillations of period 4x/(A+q) = 2nfh .. 


with slowly varying amplitudes Q cos 4(A—)t, Q sin }(4—)t. The oscilla- 
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tions of the amplitudes are 90° out of phase with one another, so that 
when one is at a maximum, the other is zero. This means that energy will 
be transferred rhythmically between the circuits. This phenomenon is a 
characteristic feature of coupled oscillatory systems and is also exhibited 
by two simple pendulums whose points of suspension are mechanically 
connected, Of course, in any actual system of this type, the circuits will 
not be entirely devoid of resistance and the oscillations will be damped 
and will ultimately die away. ia 


5.4 Alternating Currents In this section, circuits in which the applied 
¢.m.f. varies sinusoidally with the time will be studied. An e.m.f. of this type 
is supplied by an alternator, in which a coil is rotated in a magnetic field 
and its magnitude at time ¢ is given by an expression of the form 
E,, cos(wt +). E, is the amplitude, « is the angular frequency (w/2n is the 
frequency in cycles per second or Hertz) and « is the phase. It will be assumed 
that the circuit has reached a steady state in which the currents flowing 
are also oscillating sinusoidally with the same frequency as the applied 
emf. (ie, transients are neglected). 

It is convenient to regard the applied e.m.f. as the real part of the 
complex expression E, exp[j(wt+«)] = Eexp(jot), where E = E, exp(ja). 
E is termed the complex emf. Note that |E| = E, and argE = a. The 
currents are treated in the same way, a complex current I being associated 
with an actual current &I exp(jat). 

If/ is the complex current through an inductance L, the actual potential 
drop across the clement is 


LA ate) = RjoLle™. 


Thus, the complex potential drop is jwLJ. jwL is called the complex 
impedance of the element at the angular frequency . Hence, if V is the 
complex potential drop, Z is the complex impedance and / is the complex 
current, then 


V=ZI. (5.8) 


(c.f. Ohm’s Law, §3.1). A similar argument proves that equation (5.8) is 
valid for a resistance R and a capacitance C; the respective complex 
impedances are R and 1/jaC. 

A.C. circuit problems can now be solved by application of Kirchhoff’s 
two laws in the manner explained in §3.1. The currents flowing are 
represented by complex quantities and the potential drops are calculated 
from the complex impedances of the elements by use of equation (5.8). In 
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the solutions of the problems which follow, the complex currents, complex 
impedances, etc., will be referred to simply as currents, impedances, etc. 


Problem 5.6 A coil having resistance R and inductance L, and a con- 
denser having capacitance C, are connected in parallel across an alternator 
supplying an emf. Ecosat. Determine the frequency at which the 
alternator current is a minimum. 


Solution. Let 1, , 1, be the currents supplied to the coil and condenser 
respectively (Fig. 5.5), Since the impedances of these elements are R +joL, 


/ 


Fig. 5.5 


1/jaC respectively, it follows that E = (R+joL)I, = 1,/jaC. Hence, if I 
is the current supplied by the alternator, 
1—@*LC+jaCR 
R+joL 
The amplitude of the actual current is accordingly 
(1—w? LC) +.w?C?R?]* 
Ith [Sa] 


1=1,+1,=E 


and its phase is 
arg! = arg(1—w*LC+jwCR)—arg(R +joL) 


= tan-"( Cc) tan( 
1—w*LC R 
Differentiating |/|? with respect to @?, it will be found that a minimum 
exists for 


oF = ald +2a)?~a], 
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be J, Fexp(—2jn/3), Iexp(2jn/3) respectively. Since the potential drop 


== 


©» phere = CR2/I = 1/0? (Q is the quality foctor), At this frequency, the 


system is said to be in a state of amplitude resonance. across BC is £, it follows that E = ZJ. The current taken from the supply iy 


The frequency at which the current J is in phase with the applied e.m.f. at A is, therefore, 
E is determined by the equation 4 olan . 
7 oCR ob. Texp#—rexp 2 = 21 sin = ve. 
4 T=o3LC ~ RS This current has amplitude \/3E/|Z|; the currents taken at B,C have the 
ie, = (1—a)/LC. This is the state of phase resonance. If « is small, the same amplitude, but are out of phase by an angle 27/3. 3 : 
two states of resonance are indistinguishable. a In case (ii), the circuit is as shown in Fig. 5.6(ii); this is called a star 
connection of the elements Z. The currents in AO, BO,CO aretakentobe 
Problem5.7 Three-phase a.c. supply of voltage amplitude E is fed to the | J, Jexp( —2jx/3), J exp(2jx/3) respectively. By consideration of the poten- 
-yertices of a triangle. Calculate the current taken from the supply (i) when tial drops along BOC, we arrive at the equation 
these vertices are connected to one another by elements of impedance Z, 
and (ii) when these vertices are each connected to a common point O by | E = ZJlexp(—2jn/3)—exp(2jn/3)] = —/3jZJ. 
_ elements of impedance Z. Thus, J = jE/,/3Z is the current taken from the supply at A. The ampli- 
_ Solution. Three-phase a.c. supply is carried by three conductors and is tudes of the currents taken at A, B, C are all E/,/3|Z| and are mutually out 
ich that the alternating voltage across any pair of the conductors differs | of phase by an angle 2n/3. bake 
_ in phase from the voltage across any other pair by 2z/3. Thus, if A, B, C on 
are the vertices of the triangle connected to the supply, the complex Problem5.8 In the telephone circuit shown in Fig. 5.7, the inductance of — 
_ electromotive forces applied across BC, CA, AB can be taken to be E, the telephone windings is included in L’ and all resistances are negligible. 
_ Eexp(—2jn/3), E exp(2jn/3), respectively. The emf. supplied is E cos wt. Calculate the current through the telephone — 
re In case (i), the circuit is as shown in Fig. 5.6(i); this is called a delta ‘ and the frequency at which there is no signal. yi’ 
connection of the elements Z to the supply. It follows from the symmetry 5 
of the system that the currents in the arms BC, CA, AB can be assumed to 


. Fig. 5.6 Fig. 5.7 
4 


Solution. Let I be the current through the telephone and J the current 
taken from the supply. Applying Kirchhoff’ second law to the left-hand 
and right-hand circuits, we derive the equations 

E = joLJ+joMI1+(J—D/joC, 
0 = joL'I+joMJ—(J—Djjo. 
Solving for J, we obtain 
Ae (1+@?MOE 
© jo[L+L'+2M —@*C(LL’— M?)]" 

Clearly, J vanishes when »? = —1/MC; w? must, of course, be positive 
and the senses of the windings of the coils L, L’ must therefore be so 
arranged that M is negative. o 


5.5 The Lorentz Force The force F which acts upon a point charge q 
moving with velocity v in an electromagnetic field is given by 
F = q(E+vxB). (5.9) 

The definition of the electric intensity E accounts for the component gE. 
The motion of the charge creates an equivalent current element qv; 
substituting this for i ds in equation (4.6), we obtain the second component 
of force qu x B. 

F is called the Lorentz force. Note that its magnetic component is 
perpendicular to the direction of motion of g and, hence, does no work. 


Problem5.9 A point charge moves in a uniform magnetic field. Discuss 
the motion. 


Solution. Let Oxyz be stationary rectangular axes and let (x, y, z) be the 
coordinates of the charge q at time t. Take the magnetic field to be directed 
along Oz, so that B = (0,0, B). The velocity of the charge is given by 
v = (X,), 2). Hence, the force acting is given by F = q(By, — Bx,0), and, 
if m is the mass of the charge, its equations of motion are accordingly 


mx = qBy, mj = —qBX, mi =0. (5.10) 
The last of these equations shows that the motion parallel to Oz is 
uniform. The other equations can be integrated to give % = wy, = —wx, 


wherew = qB/mand the origin O is assumed to be translated appropriately 
so that the constants of integration vanish. Elimination of y now leads to 
the equation ¥+m?x = 0, and hence x = acos(wt+a). By proper choice 
of the instant t = 0, « can be made zero and then x = wy shows that 
y = —asin wt. Thus the equations of the particle’s trajectory are 

z = bt+c; 


X = acoswt, y = —asinot, 
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these are the equations of a circular helix. If the charge is constrained to 
move in a plane parallel to Oxy, b = 0 and its trajectory is a circle; if w 
is positive, the charge moves around the circle in the negative (clockwise) 
sense with angular velocity w. o 


Problem 5.10 The point charge of the last problem is constrained to 
move in the plane Oxy and, in addition to the magnetic field, is acted 
upon by an electric field of intensity E cos wt parallel to the x-axis. If the 
charge is initially at rest at O, show that its trajectory is a spiral. 


Solution. The electric field exerts an additional force (qE cos wt,0,0) on 
the particle and the equations of motion (5.10) have therefore to be amended 
to read: 

X¥ = wp+acosot, 
where @ = qB/m and « = gE/m. 

Integrating the second equation under the initial conditions x = 0, 
= 0,weget = —cxand the first equation then yields¥+w?x = acosat. 
Thisequation hasa general solutionx = Pcos wt +Q sin wt +(at sinet)/2e. 
At 1=0, x =0 and ¥ = 0; thus P = Q = 0. Substituting for x in 
J = —ox and integrating under the initial condition y = 0 at ¢ = 0, y is 
now found. The complete solution is 


yp = —@x, 


_ att 
a7) 
Putting wt = 0, «/2m* = B, parametric equations for the trajectory are 
found in the form x = fOsin 0, y = BO cosO—f sind. 

A little consideration reveals that these are the equations of a spiral. 
For large values of 0, the term f sin 0 is small by comparison with 0 cos 0 
and approximate equationsfor the trajectoryarex = B0sin0,y = BO cos 0; 
the corresponding polar equation is p= ph (p = /(x?+y?), 
@ = tan™'(y/x) = 0), which is clearly the equation of a spiral. oO 


at. a. 
x = >—sin wt, COs WI —~ ——5 sin wt. 
20 2a 


5.6 Maxwell’s Equations Thus far in this chapter, it has been assumed 
that the magnetic field due to a variable current can be calculated by the 
methods given in Chapter 4 for the case of steady currents. The implication 
is that equation (4.11) remains valid, even when j varies with ¢, Such an 
assumption is clearly unjustified since, if the divergence of both members 
of this equation is taken, the result is div j = 0 (divcurl H vanishes, 
L. Marder, Vector Fields, p. 36) and this equation is incompatible with 
equation (3.4) unless p is independent of 1 (i.e. steady conditions). If, 
however, the currents vary comparatively slowly, the term dp/ét is 
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negligible by comparison with div j, and equation (4.11) can be accepted 
as an approximation without inconsistency. We say that the conditions 
are quasi-steady. 

For general non-steady conditions, Maxwell assumed that the equations 


div D = p, div B = 0, (5.11) 
remain valid. Then, equation (3.4) can be written 
air 2) =0 (5.12) 
and this is consistent with an amended equation (4.11), namely 
,, 0D 
curl H = Ita (5.13) 


0D/ot is termed the density of the displacement current. Maxwell postulated 
equation (5.13) as a fundamental equation of the electromagnetic field. 

The set of fundamental equations was completed by the addition of one 

_ representing Faraday’s law. If C is any fixed closed curve lying in an 

~ electromagnetic field, the e.m-f. in a circuit coincident with C is $E-ds, 


ds being an element of C. Let S be any surface bounded by C and a the 
unit normal to S. Then, equation (5.1) can be expressed in the form 


d OB 
[e« = -§[ ew - [Fras 


" Applying Stokes’ theorem to the left-hand integral, we deduce the final 
Maxwell equation, an 


ih culE = —>. (5.14) 


In a solid material medium, the following equations are also usually 
valid, 


D= 8, E, B=yp,H, j=cE. (5.15) 

At a boundary between two media, assuming that the surface densities 

of charge and current are zero, Maxwell’s equations require that the 

normal components of B and D, and the tangential components of H and 
E, must be continuous. 


Problem 5.11 An electromagnetic oscillation is taking place inside a 
___ perfectly conducting spherical shell of radius a. Taking the centre of the 
___ shell as the pole of spherical polar coordinates (r, 0, @), show that Maxwell's 
equations can be satisfied by taking E, = f(r)cos 0 exp(jot), 


E, = asin Oexpljant, E, = 0, H, = Hy =0, H, = Msn Dexpljan), 


provided f, g, h are related by certain differential equations. Putting 
f =r-4y, or/e = x (c? = 1/p, &, ), Show that y satisfies Bessel’s equation _ 
and deduce that @ must satisfy the equation tan¢ = (/(1—{?), where 
¢ = wae. 
Solution, Inside the shell, j and p both vanish and B = y, H, D = &, E. 
Thus, Maxwell's equations reduce to 
div E = divH = 0, 


oH 


curl H = 9 5 curlE = —Hy ee: 


Expressing these equations in spherical polar coordinates (L. Marder, 
Vector Fields, p. 64) and substituting the assumed forms for the spheri 

polar components of H and £, it will be found that some equations are — 
satisfied identically and the remainder reduce to the forms 


d 
ar +29 =0, 
AB i 
= = joes, 


1d : 
Far h) = —jotyg, 


ld f 
par lOtT = —jop,h. 


Elimination of h between the middle pair of equations yields the first 
equation; this shows that only three of these equations are independent. 
Eliminating h from the last equation by use of the second equation, we he 
find that 

d 2.2" 

5 o+ 11-$5) =0. (5.18) 
Equations (5.17), (5.18) now give 

2 12.2" Ss 

Ser -1(2-“F) =0. (5.19) 


To solve this equation, put f = r~#y and wr/c = x; this transforms the — F 
equation to the form 


x?y”+xy' +(x? —2)y = 0; 


this is Bessel’s equation of order } (J. Heading, Ordinary Differential 
Equations, p. 82). Accepting the solution which is bounded at x = 0 
(i.e. the pole), we get 


DRL oe 
y = J,X) = () x” 3(sin x—x cos x). 


The functions, f, g, h now follow immediately: 


f= (2) g= -¥4[4(2)} h= Bory 


an arbitrary multiplier has, of course, been omitted. 

Since the shell is infinitely conducting, ¢ = « and, thus, j = cE 
implies that E = 0 within the shell. The boundary conditions over r = a 
are accordingly E, = E, = 0, i.e. g(a) = 0. This leads to the condition 


+( 2) 
=| cosx—-——] = 0 
dx x 


for x = ¢, or tan¢ = ¢/(1—27). This last equation determines possible 
frequencies of oscillation in this mode (the smallest root is € = 2:74). 
The lines of magnetic intensity are circles in planes perpendicular to 
the axis @ = 0 and the lines of electric force are curves in the meridian 
planes @ = constant. Such a device is called a cavity resonator. o 


Taking the curl of both sides of the last of Maxwell’s equations in empty 
‘space (equations (5.16)) and employing the third equation, H is eliminated 
to give curlcurlE = —y, 2, 07E/dt?. But curlcurl E = grad div E-V?E 
(L. Marder, Vector Fields, p. 39), Hence, since div E = 0, E satisfies the 
wave equation 

foe TOR 
VE= 2a (5.20) 


where c? = 1/p1)¢,. ¢ is the velocity of propagation of electromagnetic 
‘waves; substitution of the values for ¢,, 4, given in §§ 1.1, 2.5, shows that 
cis very nearly 3 x 10°ms~'. This is the velocity of light in vacuo, indicat- 
ing that light is an electromagnetic phenomenon. Elimination of E, shows 
that H satisfies the same wave equation. In uncharged, isotropic material, 
having dielectric constant ¢ and permeability x, Z and H also satisfy the 
wave equation, but the propagation velocity is a, where a? = 1/(ee, up1p). 

A plane, monochromatic (i.e. definite frequency), plane polarized wave 
being propagated in the direction of the positive z-axis with velocity a 
in a medium (e, 2) has components 


90 


E = {E, expLio(t—2/a)],0,0} 
H = {0, Hoexp[jo(t—z/a)],0}, 
where (e¢,)*E, = (ifto)*Hy (it is understood that Eo, H, may be complex 
and that real parts are to be taken). The reader should verify that Maxwell's 
equations are satisfied. A similar wave being propagated in the opposite 
sense is determined by 
E = {E, expLjo(t +z/a)], 0,0} 
H = {0, —H, expLio(t +2/a)], 0}, 
where (2¢5)#Ey = (tly)*Hy- 
Problem 5.12 The region z > 0 is occupied by a uniform conductor of 
conductivity ¢ and dielectric constant ¢ and the region z < 0 is empty. 
The wave determined by equations (5.21) is incident normally upon the 
face z = 0 of the conductor. Assuming this gives rise to a field in the 
conductor determined by 
E, = Eexp(jot—kz), E, = E, =0, 
H, =0, H,= H exp(jat —kz), H, = 0, 
and a reflected wave in the vacuum of the form (5.22) and assuming 
o > wee,, show that the ratio of the amplitudes of the reflected and 


incident waves is 
1+4- (24) 
1+A+,/(24)]° 


(5.21) 


(5.22) 


where 4 = wye,/o. 
Solution. Setting p = 0, j = cE, D = se, E, B = up, H in Maxwell's 
equations, these become 
div E = div H = 0, 
curl H = cE +26, 0E/ét, 
curl E = — ipl, 0H/At. 


These are valid in the conductor. It is found that these equations are 
satisfied by the assumed forms of E and H provided E and H satisfy 


(o+jeee,)E—kH = kE—joup,H = 0. 


For these equations to have non-zero solutions in E and H, it is necessary 
that 


aot+joce, —k 
k jot 
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aE ee ee ee ee, oa oe ee. 


ie. k? = OU (jo — wee, ). Ifo > weey, this gives k = + (Sopp, o)* (1+) 
(NB. positive root must be taken so that the field vanishes as z + + co.) 
With this value of k, the ratio H:E is determined by the equation 
HIE = (o/20yn.)}*(1—)). 
Taking the reflected wave to be given by 

E = [Ey exp{ja(t+z/c)},0,0], 

H = [0, — Hy exp{jo(t+z/c)},0], 
the net disturbance in the region z < 0 is found by superposing this wave 
‘on the incident wave (5.21) (with a = c). Thus, the field in the vacuum at 


2 = Ois given by 


E = [(E,+E,)exp(iot),0,0], HW = [0,(Hy—H')exp(jor), 0). 
The field in the conductor at z = 0 is 
E = [Eexp(jat),0,0],  H = [0, Hexp(ja), 0). 
Since E,, H, must be continuous across z = 0, the boundary conditions 
are 
E,+E,=£, H.-H, =H. 
The second of these conditions implies that 


+ 4 
fo rey af Sos ae 
(2) (Eo B= (sa) Oe 


We can now eliminate E and solve for the ratio E,:E, ; the result is 
Ey _ (2 op)*—o%(1—)) 
; E,  2ey a) +o%(1—)" 
_ The amplitudes of the incident and reflected waves are |E,| and |E’,| 
respectively; hence, taking moduli, the result stated now follows. 


In non-ferromagnetic materials, » = 1 and A will therefore be small. 
The ratio then approximates to 1—,/(2). o 


EXERCISES 


1. A current i is maintained in the long straight wire shown in Fig. 4.6 
If the centre of the circular wire is made to oscillate in the line OA with 
simple harmonic motion of small amplitude ¢ and angular frequency «, 
its mean distance from the straight wire being c, show that an alternating 
electromotive force of amplitude ji, iew[1—c//(c?—a?)] is induced 
within it. 

2. Show that the mutual inductance between a square circuit of side 2a 
and a long straight current lying in the plane of the square parallel to one 
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of its sides and distant ¢ from its centre (c > a), is 


ot og( <2) 
n c-a 


3. The primary and secondary circuits of a transformer are identical and 
have resistance R and self-inductance L, Their mutual inductance is 
L/k (k > 1). Att = 0, the primary terminals are connected to a source of 
constant e.m.f, E and the secondary terminals are short circuited. Show 
that the maximum current which subsequently flows in the secondary is 


E (k- 1 » 
R (kere: 


If the coupling is tight, so that k approximates to 1, prove that the peak 
current is almost E/2R. 


4. The primary and secondary windings of a transformer have self- 
inductances L, L’ respectively and negligible resistance. A resistance R 
and capacitance C are connected in series across the secondary terminals 
and an emf. of amplitude E and angular frequency q is applied to the 
primary terminals. If M is the mutual inductance between the windings, 
show that the current in the secondary is a maximum when w* = 
L/C(LE—M?) and that the equivalent impedance across the supply is 
then LR'/M. 
5. The space between two very long, circular, coaxial cylinders is 
evacuated and electrons are set free with negligible initial velocity at the 
surface of the inner cylinder. They are accelerated from the inner cylinder, 
of radius a, towards the outer, of radius b, by means of a constant potential 
difference V between the cylinders. Ifa uniform magnetic field B is applied 
parallel to the axis of the cylinders, show that the electrons will just reach 
the outer cylinder if ‘ 
8mb?V = eB(b? —a?)?, 
e being the charge and m the mass of the electron. (Magnetron) 


6. Repeat the cavity resonator calculation of Problem 5.11 for an oscilla- 
tory mode in which 

H, = f(rcosOexp(jot), Hy = g(r)sin@exp(jmt), Hy = 9, 

E, = E, = 0, Ey = h(r)sin 0 exp(ja) 


and show that the possible frequencies @ are determined by the equation 
tanf=¢. 
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Appendix 


International System of Units 
(Rationalized MKSA System) 


Quantity ‘SI unit 
- Length metre (m) 
Mass kilogramme (kg) 
Time second (s) 
Force newton (N) 
Work and energy joule (J) 
Power watt (W) 
coulomb (C) 
ampére (A) 
I volt (V) 
_ Electric dipole moment coulomb-metre 
Electric displacement D coulomb-metre~? 
_ Electric intensity E newton-coulomb~? or yolt-metre~! 
" Electric polarization P coulomb-metre~ 
Flux of magnetic induction N weber (Wb) 
Magnetic induction B weber-metre~? or tesla (T) 
a Magnetic intensity H ampére-metre~! 
a Magnetic pole strength p ampére-metre 
_ Magnetic moment m ampére-metre? 
Magnetization M ampére-metre~* 
fs _ Capacitance € farad (F) 
Resistance R ohm (Q) 
Inductance L henry (H) 
Frequency v = w/2n hertz (Hz) 


| | 


Accumulation of flowing charge 58 
Alternator 82 

Alternating emf. 76, 82 

Ampére’s Law 72, 73, 74 
Amplitude 82 

Angular frequency 82 

Anode, cathode ray 7 


Back em. 76 

Bessel's equation 90 

Biot-Savart law 63 

Boundary conditions, current flow 57 
Rome conditions, electrostatic 12, 


Boundary conditions, electromagnetic 
88 
Boundary conditions, magnetostatic 43 


Capacitor 25, 79, 82 
Capacitance 26 

Cavity resonator 90, 94 

Charge inside earthed hemisphere 23 
Charge near earthed 20 

Charge near uncharged sphere 21 
Circular current, field of 63 
Condenser 25 

Conductor, charge on surface of 12, 32 
Conductor, no charge in interior of 10 
Conductor, Posen of 10 
Conductivity 54 

Conservative field 3 

Contact potential $7 

Continuity, equation of 54 

Complex current 82 

Complex emf. 82 

Coulomb's law 1, 21, 31 

Couple on plane current loop 69 
Coupled oscillatory circuits 80 
Current density 54 

Current density, flux of $4 

Current loop as dipole 70 

Current strength S51 

Cylinder, magnetized 44 

Cylinder near to earthed plane 26 
Cylindrical capacitor 26, 28 

Cylinders, electric field between two 25 
Cylinder, field of charged 9 


Dampled oscillations 82 


Datum for potential 4 
Delta connection 84 

Demagnetizing field 45 

Dielectric 29 

Dielectric between capacitor plates 32 
Dielectric constant 30 

Dielectric cylinder and line charge 38 
Dielectric plane slab and point charge 36 
Dielectric sphere in uniform field 34 
Disc, uniformly charged 1 
Displacement current 88 
Displacement, electric 31 


Dipoles, forces between 45 
Dipole in spherical cavity 16, 35 
Dipole near carthed sphere 22 
Divergence theorem 11, 54 
Doublet, line 27 

Doublet, moment of 27 


Earth, potential of 10 
Earthed conductor 10 
Electrode 5S 

Electrode, circular 55 
Electrode, cylindrical 55 
Electrode, spherical 55 
Electromagnetic field, steady 63 
Electromagnetic waves 90 
Electromotive force 51 
Elliptical current, field of 75 
Energy of capacitor 26 
Energy of magnetic dipole 45 
Equipotential 12 


Faraday's Law 76, 79, 88 

Ferromagnetic materials 42 

Flux of electric intensity 8 

Flux of induction 69, 76 

Force between straight currents 67, 71 

Force between straight and circular 
currents 68, 75 

Force between plates of capacitor 17 

Force on conducting surface 17 

Force on current 67 

Frequency 82 


Gauss’s law 8, 12, 21, 26, 31, 32, 41 
Gradient of potential 4 


oe 


Velocity of light 90 
Waves, plane 90 
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